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ABSTRACT 

We have carried out a survey for CO J=l-0 and J=2-l emission in the 260 early-type galax- 
ies of the volume-limited ATLAS 3D sample, with the goal of connecting their star formation 
and assembly histories to their cold gas content. This is the largest volume-limited CO sur- 
vey of its kind and is the first to include many Virgo Cluster members. Sample members are 
dynamically hot galaxies with a median stellar mass ~ 3 x 10 10 Mq; they are selected by 
morphology rather than colour, and the bulk of them lie on the red sequence. The overall CO 
detection rate is 56/259 = 0.22 ±0.03, with no dependence on K luminosity and only a mod- 
est dependence on dynamical mass. There are a dozen CO detections among the Virgo Cluster 
members; statistical analysis of their H2 mass distributions and their dynamical status within 
the cluster shows that the cluster's influence on their molecular masses is subtle at best, even 
though (unlike spirals) they seem to be virialized within the cluster. We suggest that the cluster 
members have retained their molecular gas through several Gyr residences in the cluster. There 
are also a few extremely CO-rich early-type galaxies with H2 masses > 10 9 Mq, and these 
are in low density environments. We do find a significant trend between molecular content and 
the stellar specific angular momentum. The galaxies of low angular momentum also have low 
CO detection rates, suggesting that their formation processes were more effective at destroying 
molecular gas or preventing its re-accretion. We speculate on the implications of these data for 
the formation of various sub-classes of early-type galaxies. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: ISM — 
galaxies: structure — Radio lines: galaxies. 
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1 INTRODUCTION 



Understanding galaxy formation and evolution is at the heart of 
much of current astrophysics. The task is, however, greatly hin- 
dered by the lack of quantitative, physically-driven (as opposed to 
empirical) models of the regulation of molecular gas and star for- 
mation on galactic scales. All galaxies or protogalaxies must have 
begun as gas-rich entities, and during their first epochs of star for- 
mation would have belonged to the "blue cloud" in a galaxy colour- 
magnitude diagram. In contrast, today there is a clear separation 
between the blue cloud (composed mainly of star-forming disc and 
dwarf galaxies) and the red sequence (old and dynamically hot stel- 
lar systems). Thus, an outstanding question in galaxy evolution is 
how the present-day early-type galaxies moved quickly from the 
blue cloud to the red sequence. 

This rapid movement to the red sequence requires, at mini- 
mum, an abrupt cessation of star formation at high redshift so that 
the g lobal colours change from blue to red (IThomas et al . 2005, 
2010). Since today's early-ty pe galaxies are generally rather poor in 
atomic and molecular gas dLees et al.ll99ll) , it is often assumed that 
the cessation of star formation was achieved by removing, destroy- 
ing, or consuming the cold gas. The most massive early-type galax- 
ies are often rich in hot gas (lO'Sullivan et al.l200ll) . so it is possible 
that molecular gas is destroyed by heating. An intriguing alternative 
suggestion is that the cold gas might be retained but re ndered un- 
suitable for star formation activity dMartig et al.ll2009l) . Thus, the 
story of the development of the red sequence is, to a large extent, a 
story of what happens to the cold gas in galaxies. 

Several mechanisms have been proposed to deplete the cold 
gas from galaxies and move them onto the red sequence. For exam- 
ple, a major merger between two gas-rich galaxies could (depend- 
ing on geometry) result in a large-scale loss of angul ar momentum , 
dropping gas to the center of the merger remnant (Barnes l Eool . 
The gas could then be funneled into a black hole and/or consumed 
in a LIRG-type burst of star formation activity. Alternatively, the en- 
ergy input into the interstellar medium from an active nucleus or a 
starburst could hypothetically destroy and/or unbind the cold gas, as 
well as prevent hot gas from cooling. Environmental effects such as 
ram pressure stripping, gravitational interactions and galaxy harass- 
ment could also be important in at least some early-type galaxies. 
In fact, it is likely that all of these processes are important at some 
level, but we do not yet have strong constraints on their relative im- 
pacts. The behavior of the cold gas in all of these transformational 
processes needs to be better understood. 

It has also been known for some years that today's early- 
type galaxies, while relatively poor in cold gas compared to 
spirals, are not completely devoid o f cold gas. The gas is de- 
tected in dust y silhouette discs (e.g . iGoudfrooii et al.lll994l). HI 
emission (e.g.lWardle & Knaprj[l 986; Huchtmeier, Sage & Henkel 



ll995t lMorganti et al.ll200d;IOosterloo et al.ll2010h. and CO emis 



IWiklind et al.lll995|- ICombes etalj 120071 : 
Welch, Sage, & Young |2010|) . The origin of this gas is not well 
unde rstood; it could have come from internal stellar mass los s 
(e.g. iMafhews & Brighentil 120031 : ICiotti. Ostriker & Progal l20ld) , 
or it could have been acquired from an external source such 
as another galaxy or cold mode accretion. In some early-type 
galaxies the case for an external origin is very clear because the 
specific angular momentum of the gas is dram atically different 
from that of the stars. Cen A dOuillen et al.lll992l) and NGC3032 
( Youn g. Bureau & Cappeilari| [2008) exhibit this property, but in 
many other cases the gas kinematics are consistent with an internal 
origin. Again, probably both external and internal processes con- 



sion dLees et al. 1991 
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tribute at some level. A quantitative understanding of their contri- 
butions would make it possible for the gas to serve as an indicator 
of the galaxies' evolutionary histories and the processes that moved 
these galaxies onto the red sequence. 

These considerations motivate a study of correlations between 
the cold gas content of early-type galaxies and their other proper- 
ties. For example, the kinematic structures in early-type galaxies 
(dynamically cold stellar discs, counterrotating cores, or even the 
absence of a measurable net rotation) preserve clues to their assem- 
bly histories. If the cold gas in these galaxies is leftover from the 
time of their assembly, we might expect the gas content to be cor- 
related with kinematic structure. In addition, it would be useful to 
test the influence of AGN activity and/or local environment on the 
gas content. Looking to the future, observations of the cold gas in 
early-type galaxies can also be used for testing theoretical models 
of the ISM structure and of star formation processes. 

We address these goals through a CO survey of the ATLAS 3D 
sample, a large and complete volume-limited sample of early-type 
galaxies spanning a variety of environments and AGN activity lev- 
els. This is the largest survey to date of CO in early-type galaxies, 
being factors of 5 to 10 larger than previous studies with similar 
sensitivity. It is also the first volume-limited sample to have spa- 
tially resolved stellar kinematics, stellar populations, and ionized 
gas data for all target members. And while the smaller C O survey 
of the SAURON early-type galaxies dCombes et alj|2007h included 
both Virgo Cluster members and field galaxies, it was not a volume- 
limited sample. The present paper is the first to allow reliable statis- 
tical analyses of gas contents inside and outside the Virgo Cluster. 

Section [2] summarizes the selection criteria for the ATLAS 3D 
sample. Section[3]describes the 204 new CO observations obtained 
for this survey and the literature data collected for the remainder. 
Section [4] describes global results of the survey, and section [5] an- 
alyzes the molecular gas content compared to a variety of galaxy 
properties including stellar kinematics and environment. Correla- 
tions with stellar populations, star formation activity, and AGN ac- 
tivity will be discussed in a future paper. The discussion in section[6] 
covers implications for the origin of the molecular gas in early-type 
galaxies, the removal of gas in clusters and the possible transforma- 
tions of spirals into early-type galaxies. 



2 SAMPLE 

The ATLAS 3D sample is a complete volume-limited sample of 
early-type galaxies brighter than Mk = —21.5, covering distances 
out to 42 M pc, with some restriction s on Declination and Galactic 
latitude (see Cappella ri et al.ll2011al hereafter Paper I). The early- 
type sample is actually drawn from a parent sample which has no 
colour or morphological selection, and optical images of the en- 
tire parent sample have been inspected by eye for large-scale spiral 
structure. The presence or absence of spiral structure is, of course, 
a physically motivated selection criterion because it tests whether 
the stellar population is dynamically cold (and self-gravitating) or 
not. The 260 galaxies lacking spiral structure form the basis of the 
ATLAS 3D project. Integral-field optical spectroscopy was obtained 
with the SAURON instrument on the William Herschel Telescope; 
it covers a field of at least 33" x 41", typically extending to one 
effective (half-light) radius R e . 

Paper I gives additional information on the details of the mor- 
phological selection of the sample, assessment of Virgo Cluster 
membership, the distances of the galaxies and measured optical ve- 
locities. Paper I shows that the morphologically-selected ATLAS 3D 



galaxies clearly trace out the red sequence, with a small population 
of bluer galaxies in th e "transition region" be tween the blue cloud 
and the red sequence. iKrainovic et al.l d201 lh (Paper II) present a 
morphological classifica tion of the internal stel lar kinematic struc- 
ture of the galaxies, and lEmsellem et all J201 lh (Paper III) analyse 
the stellar specific angu lar momentum proxy Xr for the galaxies. 
Cappe ilari et al.l l l2011bh (Paper VII) present measurements of the 
local galaxy density, which we use for study of environmental ef- 
fects on cold gas content. 



3 DATA 

3.1 IRAM 30 m observations and data reduction 

The IRAM 30 m telescope at Pico Veleta, Spain, was used for si- 
multaneous observations of 12 CO (1-0) and (2-1) during July 2007, 
March 2008, and November 2008. The beam FWHM is respec- 
tively 21.6" and 10.8" at the two frequencies. The SIS receivers 
were used for observations in the wobbler switching mode, with 
reference positions offset by ±100" in azimuth. The 1 MHz fil- 
terbank back-end gave an effective total bandwidth of 512 MHz 
(1330 km s _1 ) and a raw spectral resolution of 2.6 km s 1 for 
CO(1-0). The 4 MHz filterbank gave an effective total bandwidth 
of 1024 MHz (also 1330 km s -1 ) and a raw spectral resolution of 

5.2 km s' 1 forCO(2-l). 

If the CO emission is quite extended, the relatively small wob- 
bler throw could potentially result in some emission being located 
in the "off' positions. The signature of this problem would be a neg- 
ative feature, most likely at a velocity different to the optical sys- 
temic velocity. There is no strong evidence of this problem in any 
spectra except fo r those of NGC 464 9, where known CO emission 
from NGC 4647 dYoung et al.ll2006l) appears as a negative feature 
offset by 270 km s _1 . The effect may possibly be seen in the spec- 
tra of NGC 2685, which however is already known to have extended 
CO emission (see below). 

The system temperatures ranged between 190 K and 420 K at 
2.6 mm, and between 240 K and 600 K at 1.3 mm. A small number 
of targets had significantly poorer quality data, with system tem- 
peratures up to 1400 K for CO(2-l), and for NGC 2698 the system 
temperatures were so large at 2-1 as to make those data useless. 
The pointing was checked every 2 to 3 hours on a nearby planet 
or bright quasar, and focus was checked at the beginning of each 
night as well as after sunrise or more often if a suitable planet was 
available. The time on source typically ranged from 12 minutes to 
36 minutes and occasionally longer, being weather dependent, in- 
teractively adjusted so that the final co-added 12 CO (1-0) spectrum 
for each galaxy had a rms noise level near 3.0 mK (TX) ~ 18 mJy 
after binning to 31 km s _1 channels. 

The individual 6 minute scans were all inspected and those 
with unusually poor quality baselines or other problems were dis- 
carded, to reduce the possibility that a faint emission line feature 
could be created by baseline instabilities. The good scans were aver- 
aged together, weighted by the square of the system temperature. In 
six cases the CO(2-l) spectrum required a linear (first order) base- 
line, while in all other cases a simple constant (zero order) baseline 
was used. If no emission was apparent in the spectrum, all of the 
channels were used to determine the baseline level. If emission was 
obvious the baseline was determined using interactively selected 
line-free regions of the spectrum. 

Integrated intensities for each galaxy were computed by sum- 
ming the spectrum over velocity. In the cases with weak or absent 
line emission, the velocity range is 300 km s _1 centered on the 
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systemic velocit y of stellar absorption lines as measured from our 
S AURON data ( Cappellari et alfeoilah . In most cases our new ve- 
locities agree with values listed in the HyperLEDA database. For 
NGC4486A our new velocity disagrees by 700 km s" 1 due to the 
presence of a bright foreground star corrupting previous long-slit 
spectroscopy, and our original use of the LEDA velocity means 
that we have effectively no CO data for that galaxy. For obvious 
detections the integrated velocity range covers the real emission. 
Table [4] lists the rms noise levels in each spectrum, the velocity 
ra nges summed, a nd the integrated intensities h-o and I2-L As 
in ISageetaljj2007h . the statistical uncertainty 07 in a sum over JVj 
channels, each of width Aw and rms noise level a, is 



aj = (A«)VjV,(l + JV»/JVj). 



(1) 



A contribution from the uncertainty in estimating the baseline level 
is included (Ni /Nb is the ratio of the number of channels in the line 
to the number of channels used for measuring the baseline level). 
When the integrated intensity is greater than three times its own 
uncertainty the galaxy is counted as a detection. Detected lines were 
also fit with Gaussian or double-horned template profiles. The fitted 
line areas are consistent with the summed areas; specifically, for 30 
galaxies with 1-0 integrated intensities > 3 K km s _1 , the mean 
ratio (sum area)/(fit area) is 0.96 with a dispersion of 0.04. The 
properties of the line fits are discussed in greater detail in section 
|4~2lbelow. 

Figure[T]shows the integrated intensities of all sample galaxies, 
in units of their own statistical uncertainty (equation [TJ. The long 
positive tail represents real detections; the peak symmetric about 
0.0 is caused by thermal noise and low level baseline wiggles in the 
spectra. In general, one would expect the statistical uncertainty of 
equationQjto underestimate the true uncertainty in a sum because it 
does not account for correlated channels (baseline ripples), but Fig- 
ure Q] suggests this effect is small since the dispersion in the inte- 
grated areas is not substantially larger than the theoretical estimate. 
In our sample there are a total of 56 galaxies with line intensities 
above the 3a level; four of them have 3.0 ^ I1-0/&1 ^ 4.0 and 
might be treated with caution. Of these four, NGC 4684 is also de- 
tected independently in the CO(2-l) line at a higher level of signifi- 
cance. In 203 empty spectra we would expect to find only 0.26 false 
positives with J/aj > 3.0 and also 0.26 false negative "detections" 
(noise that looks like an absorption line) with J/cj < —3.0. As 
usual, the real outlier detection rates are probably larger than that, 
but the fact that we find so few 3(J negative "detections" (Figure[T) 
suggests that there are also few false positives. 

Conversion to main beam brightness temperatures is achieved 
by dividing the antenna temperatures by the ratio of the beam 
and forward efficiencies, B c fi/F c g = 0.78 at 2.6 mm and 0.63 at 
1.3 mm. Temperatures quoted in the tables are main beam bright- 
ness temperatures (T m b). The conversion from main beam tempera- 
ture to flux density in ly is 4.73 Jy K~ 1 for the 30m at both frequen- 
cies. The total mass of molecular hydrogen is estimated using a ratio 
N(H 2 )/Ii-o = 3 x 10 20 cm" 2 (Kkms -1 ) -1 toickman et all 19861 ; 
IStrong et all 1988ll2004h . This ratio converts an integrated intensity 
to an H2 column density averaged over the beam, and multiplication 
by the beam area gives the total H2 mass as 



M^ = 6.0X10 4 
M yMpc 



J T mb dv 
K km s" 1 




A beam size of 21.8 is assumed and the effective beam area is 
1.13 9p WHM . The molecular mass thus only refers to the mass 
within the central 22" of the galaxy, which corresponds to a di- 
ameter of 2.2 kpc at a typical distance of 20 Mpc. The masses will 



2 4 6 8 10 

l 1-o/ (7 1-0 

Figure 1. Integrated intensities. For each target, the integral of the CO(1-0) 
line over the relevant spectral window (300 km s - 1 or the detected emis- 
sion) is plotted in units of its own statistical uncertainty. Detected galaxies 
have a long tail up to 7i_o/o"i-o = 42, outside the range of this figure. A 
Gaussian with a dispersion of 7i_o/°"l-0 = 1-09 is overlaid. 



be underestimated if the molecular gas is more extended than the 
beam of the 30m telescope. Based on previous experience with in- 
terferometric imaging of CO in early-type galaxies, t his effect can 
occasionally be as large as a factor of two ormore (e.g. lYoung 2002, 
l2005l ; lYoungT Bureau & Ca ppellari 2008) but it cannot be predicted 
a priori. The beam of the 30m telescope for CO(1-0) does cover 
roughly one third of the SAURON field-of-view (33" x 41" in low- 
resolution mode), so in most cases the mean stellar parameters are 
determined over a region of similar size to the CO(1-0) beam. Inter- 
ferometric observations of the CO detections are being made with 
the CARMA and Plateau de Bure arrays, and these will find molec- 
ular gas which is more extended than the 30m beam. 



3.2 Duplicate observations 

NGC 3412, NGC 3941, NGC 4026, and NGC 4346 were observed 
for this project and w ere also observed with the 30m telescope by 
IWelch & Sagel (2003). The two sets of observations have similar 
noise levels and the results in all four cases are consistent between 
the two papers. NGC 3941, NGC 4026, and NGC 4346 are non- 
detect ions in both papers. In the case of NGC 34 1 2, IWelch & Sagel 
j2003l) claim a detection in CO(1-0) at about the 3.8<r significance 
level. We find no detection but have somewhat worse noise so that 
the flux upper limit we claim is w 50% higher than the flux they 
quote. Similarly, NG C 7332 and NGC 7457 were also observed by 
IWelch & Sage] 120031) although using the NRAO 12m telescope, and 
the non-detections reported here are consistent with those data (our 
sensitivity is worse for NGC 7457). 

NGC 3193, NGC 4494, and NGC 469 7 have CO spectra from 
the 30m telescope both in this paper and in ISageetalJ0OO7h . with 
similar noise levels. NGC 821, NGC 3226, NGC 3377, NGC 3379, 
NGC 3640, NGC 4697, a nd NGC 5845 have CO s pectra published 
both in this paper and in IWelch. Sage. & Yound bOlfj) . Seven of 
these 10 (all except NGC 3640, NGC 4494 an d NGC 3226) are 
non-detections in both datasets. For NGC 4494, ISage et al.l J2007h 
claim a detection at about the 3a level in CO(1-0); the fea- 
ture is weak and broad, and is not duplicated in their 2-1 spec- 
trum. We do not find any suggestive features in our spectrum and 
so we claim a non-detection at a flux level ab out 30% smaller 
than the detection claimed by ISage et al.l fe007l) . For NGC 3640, 
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IWelch. Sage. & Yound d2010h find a 3a detection in CO(2-l) but 
not in CO(l-O); we find no emission in somewha t more sensitive 
data. For NGC 3226, IWelch. Sage. & Yound d2010l) again claim a 
3a detection in CO(l-O) at a velocity of 1100 km s _1 ; however, 
our new optical spectr oscopy indicate s the s tellar systemic veloc- 
ity to be 1321 km s~ 1 ICappellari et aL I d2011ah and our formal sum 
over 300 km s -1 centered on the optical velocity finds no detection. 
Given the nature of 3a detections, there is no glaring inconsistency 
between these observ ations. 

As mentioned in lCombes et al.l d2007» . there are a small num- 
ber of detections which are observed by tw o teams. NGC 4150 is 
detected both bv lCombes et"al] J2007h and bv lWelch & Sagel (2003) 
using the 30m; the CO(1-0) fluxes are entirely consistent, but the 
CO(2-l) fluxes differ by a factor of two. The diffe rence could easily 
arise f rom the telescope pointing, as Young, Bur eau & Cap pellari 
(2008) have shown that the CO emission in NGC 4150 is compact 
and a small pointing offset could mean that much of it is observed 
at a lower gain. NGC 3489 shows a factor of two difference in both 
CO lines between those papers, which could again be a pointing 
problem or inaccurate spectral baselines. Finally, it is worth remem- 
bering that the absolute calibration of these millimetre observations 
is not normally better than 10%-20%. These comparisons suggest 
that in the fainter detections the flux levels (hence H2 masses) might 
really be uncertain by up to 50%. This is especially true for the 
CO(2-l) line where the telescope pointing becomes both more crit- 
ical and more difficult because of the smaller beam. 

Eleven of the ATLAS 3D members have also been observed 
with the BIMA, CARMA, and Plateau de Bure interferome- 
ters. Comparisons of the single dish and interferometer total 
flux measurements for NG C 3032, NGC 4150, NGC 4459 an d 
NGC 4526 are prese nted by lYoung Bureau & Cappellaril d2008l) ; 



iung, Bureai 

for NGC 2768 b y ICrocker et all ilOofir . and for NGC 4550 
by ICrocker et al. I d2009al) . Recent observations of NGC 524, 
NGC 3489, NGC 4477, an d NGC 7457 at the Plateau de Bure are 
reported by ICrocker et all (l2010h . In general there is reasonable 
agreement, given the cautionary note mentioned above. Our 30m 
nondetection of NGC 7457 is confirmed at the Plateau de Bure. 
NGC 4550, which is a non-detection here, does have a small amount 
of molecular gas that is quite compact and more easily detected by 
an interferometer with a smaller synthesized beam. The interfero- 
metric flux of NGC 3489 is larger than what was found with the 
30m because the molecular gas is extended over an area somewhat 
larger than the 30m beam. Other detections are consistent in total 
flux and linewidth between the single dish and the interferometer. 

3.3 Literature Data 

NGC3607, NGC4111, NGC4143, NGC 4203, NGC 4251, and 
NGC 5 866 were observed with the 30m telescope bv lWelch & Sagel 
d2003l) . We use molecular masses derived from their reported CO( 1 - 
0) line intensities, scaled to the con version factor and the distance 
that we assume. The method used bv lWelch & Sagel d2003h for esti- 
mating uncertainties and upper limits is similar to the method used 
here, with the slight difference that the velocity range used for sum- 
ming varies between 200 and 400 km s _1 instead of our constant 
300 km s" 1 . Similarly, N GC 3073, NGC 36 05, NGC 4283, and 
NGC 4636 were observed bv lSageetal.ld2007l) and we take results 
from that paper. The CO emission in NGC 5866 and NGC 3607 
seems to be extended, as the line is detected in pointings offset from 
the galaxy nucleus. However, for consistency with the rest of our 
observations, only the data from th e central pointings are u s ed her e. 
Data for NGC 3522 are taken from lWelch. Sage. & Young] d2010l) . 



The most recent single dis h CO spectrum of NGC 4476 is pub- 
lished bv lWildindet1ii]|l995h . but the signal-to-noise ratio in that 
spectrum is q uite low. We consider the CO(1-0) flux measured by 
lYound 120021) in a map made with the BIMA array to be more reli- 
able and have used that value. 

NGC 4278 was obse rved at the IRAM 30m telesco pe by both 
ICombes et ail d2007h and lWelch. Sage. & Yound d201Ch . The first 
paper reports a tentative detection of 2.1±0.41 K km s _1 in the 1-0 
line, whereas the second reports a 3a upper limit of 1.9 K km s _1 . 
An observation with the IRAM P lateau de Bure inter ferometer re- 
ports a m ore sensitive upper limit ( Crocker et al The nonde- 
tection of lWelch. Sage. & Yound OOld) is used here. 



4 RESULTS 

4.1 CO detection rate 

The derived H2 masses and limits for all ATLAS 3D sample mem- 
bers are listed in the last column of Table|4] For non-detected galax- 
ies, the quoted mass limits correspond to three times the statisti- 
cal uncertainty in an integral over the assumed line width. There 
are 56 detections out of 259 sample members with CO data, giv- 
ing an overall detection rate of 0.22 ± 0.03. The quoted error is 
a la formal uncertainty assuming a binomial distribution. Thus, 
the current detection rate i s cons istent with the value 0.28 ±0.08 
reported bv lCombes et al.l d2007t) for the much smaller and some- 
what closer SAURON sample. Of the 56 total detections, 46 are de- 
tected at greater than 3a in both the CO(1-0) and CO(2-l) lines. Six 
additional galaxies (NGC 0509, NGC 2685, NGC 3245, NGC 4283, 
NGC 4476, and NGC 6798) are detected in only the CO(1-0) line, 
though NGC 4476 does not have high quality data in the 2-1 line. 
Three galaxies (NGC 3599, NGC 4036, and NGC 4643) are for- 
mally detected only in the CO(2-l) line. 

The detection rate reported here is a strict lower limit to the 
true incidence of molecular gas in early-type galaxies. The case of 
NGC 4550 was already mentioned above. In addition, NGC 2685 
(detected here in the 1-0 line but not 2-1) is known to have a sig- 
nificant amount of mo lecular gas that is located ou tside of our cen- 
tral pointing position dSchinnerer & Scovillel2002l) . For the sake of 
consistency and homogeneity, we use our own data for these two 
galaxies in the current analysis and simply bear in mind that we 
have not yet found all of the molecular gas in early-type galaxies. 
Due to the finite sensitivity of the observations there is, of course, 
also a detection bias against broad and weak lines. 

Our current detection rate is stri kingly inconsis t ent wi th the 
0.78 ± 0.17 detection rate quoted by IWelch & Sagel J2003h for a 
volume-limited sample of SO galaxies. Even if we exclude from our 
ATLAS 3D sample the 25% with types T < -3.5 (ellipticals, as 
classified in LEDA), our detectio n rate is still only 0.28 ± 0.03. The 
sample o f IWelch & Sagel d2003T) has a factor of 10 fewer galaxies 
than are presented here, but t he reasons for the dif ference in detec- 
tion rate are still not obvious. IWelch & Sagel d2003l) exclude known 
Virgo Cluster members, though as we will show later, our detec- 
tion rate is not much lower in the Virgo Cluster. We will also show 
that our detection rate is not a strong function of stellar luminosity 
over the range covered here (—26 < Mk < —21.5). Thus, other 
factors such as selection criteria and sensit ivities may contribute to 
the difference. The luminosity selection of IWelch & Sagel d2003h is 
done on a B luminosity, which means that at the faint end of the 
sample the members will be preferentially bluer than other galax- 
ies of similar stellar mass. They could thus have larger molecular 
content per unit stellar mass. And as we have mentioned above, it is 



© 0000 RAS, MNRAS 000, 000-000 



6 L. M. Young et al. 



clear that greater sens itivity would produce a greater detection rate. 
IWelch & Sae3 J 20031) have integrated down to a mass limit which 
is a fixed fraction of the optical luminosity of the galaxy, and in 
some cases this limit is several times deeper than we have been able 
to achieve. On the other hand, our limits are deeper for other ob- 
jects especially at high luminosities. Since the integra ted CO prop- 
erties are not trivially related to the optical properties (ICombes et alJ 
2007), it is not clear which of the two strategies (fixed thermal noise 
level or fixed CO/La' ratio) would be expected to produce a greater 
detection rate. 

4.2 Distribution and temperature of molecular gas 

Figure[2]shows the 30m spectra obtained in the CO(1-0) and CO(2- 
1) lines towards the detected galaxies, along with their optical ve- 
locities. The detected lines were fit with a Gaussian to provide esti- 
mates of the peak intensity, central velocity, and FWHM (Table [5). 
In a few cases the detected line is obviously better described by a 
double-horned spectrum than by a Gaussian, and these were fit with 
a template which uses a parabolic shape between two sharp line 
edges. Comparisons of Gaussian and double-horned fits revealed 
that the two shapes gave almost identical estimates for the CO sys- 
temic velocity and line width, so those parameters are robust to the 
details of the line shape in these cases of moderate signal-to-noise 
ratio. In many other cases the line is too asymmetric to be well fit 
by either of these parametrizations, and in these cases the peak in- 
tensity, central velocity, and FWHM are measured directly from the 
spectra. 

Galaxies that were fit with a double-horned profile shape in- 
clude NGC2764, NGC3182, NGC3665, NGC4429, NGC4435, 
NGC4526, and UGC 06176. Other galaxies that show signs 
of a double-horned shape such as a flat top and shoulders 
above the Gaussian, or steeper sides than the Gaussian, include 
IC 1024, NGC 524, NGC 3626, NGC 41 19, NGC 4281, NGC 4596, 
NGC4643, NGC 6014, NGC 7465, PGC058114, and UGC 05408. 
The classic mechanism for producing a double-homed shape is that 
in a disc with a flat rotation curve, the material along the major axis 
(near the line of nodes) all ha s much the same pro jected velocity 
and it is the extreme velocity dWiklind et al.lll997l) . Thus we find 
that almost a third of the detections show signs that the gas is in a 
regular disc with a flat rotation curve. Profiles that do not show a 
double-horned shape could indicate that the gas does not extend as 
far as the flat part of the rotation curve, that the disc is face-on, the 
gas is not in a disc, or that the signal-to-noise ratio of the line is 
simply too low. Thus, the true incidence of regular discs is likely to 
be larger than 30%. 

At least 13 of the detected lines are visibly asymmetric, 
and again that must be a lower limit since the detectability 
of asymmetries requires relatively high signal-to-noise ratios. 
These asymmetries could arise from a gas distribution which is 
intrinsically lopsided (as if the gas were recently acquired and 
has not yet settled into dynamic equilibrium). They could also 
arise if the gas's spatial extent is comparable to or larger than the 
size of the beam and pointing errors cause some of the gas, even 
in a regular disc, to fall outside of the beam. Notable examples 
of asymmetric profiles include IC 1024, NGC 1222, NGC 3626, 
NGC 3665, NGC 4335, NGC 4573, NGC 4694, NGC 47 10, 
NGC 5173, NGC 7465, PGC 056772, PGC 058114, UGC 05408 , 
and UGC 09519. NGC 3607 and NGC 5866 JWelch & Sagell2003l) 
also fall into this category, and their gas is known to be spatially 
extended from detections in multiple pointings. 

Systemic velocites for the CO lines (Table |5j are in good 




1 10 
l,_ (K km s" 1 ) 



Figure 3. Integrated intensities of CO(1-0) and CO(2-l) for the detected 
galaxies. For those detected in one transition but not the other, an arrow 
is plotted with the tail at the 3<r formal limit in the non-detected transition. 
Two indicative lines are also plotted, the bottom one showing 7i_o = I2—1 
and the top showing I2— 1 = 47i_q. 



agreement with our newly measured optical velocities; the dif- 
ference between those two velocities has a dispersion of only 16 
km s _1 . Outliers, for which the velocity difference is greater than 
50 km s _1 , include only cases in which our central pointing is not 
expected to have recovered all of the CO emission (NGC 2685) 
or cases in which the CO line is faint, so that it might be intrin- 
sically asymmetric or double-horned (NGC 6798, NGC 3245). If 
there were members of the sample that had recently acquired their 
gas from an external source, they might have revealed themselves 
via molecular gas which had not yet settled into the center of the 
galaxy's potential. However, we find no convincing evidence for 
any of these cases. More detailed probes of the dynamical status of 
the molecular gas will require interferometric CO maps. 

The integrated CO(2-1)/CO(1-0) emission ratio is listed in Ta- 
ble [4] and Figure [3] plots the integrated CO(2-l) intensity against 
the CO(1-0) intensity when at least one line is detected. Since only 
one central beam has been observed, the spatial region covered is 
approximately four times smaller in CO(2-l) than in CO(1-0) and 
the line ratio is affected by both the excitation temperature and the 
spatial distribution of the gas. If the CO emission uniformly fills 
both the CO(1-0) and the CO(2-l) beam and both lines are opti- 
cally thick, with the same excitation temperature, one expects the 
integrated intensity in K km/s to be the same in both lines. This re- 
sult occurs because, by definition, the measured brightness temper- 
ature is a specific intensity averaged over the beam area. If, how- 
ever, the CO emission is compact compared to the beams (point 
spread functions) then the measured intensity in the CO(2-l) line 
should be larger by up t o a factor of 4. Subthermal excitation (e.g. 
Braine & Combes 1992) would decrease the intensity of the J=2-l 
line relative to J=l-0, if densities are not large enough to populate 
the upper rotational levels. 

Most data points in Figure[3]have I2-1/I1-0 between 1 and 4, 
suggesting that in most cases the emission does not fill both beams 
uniformly nor is it extremely compact. Intermediate size scales, in 
which the extent of the emission is comparable to or smaller than 
the beams, are probably the most common. Two interesting out- 
liers are NGC 1266 and NGC 4684, which have quite high CO(2- 
1)/CO(1-0) line ratios 3), suggesting compact gas distributions. 
A few cases show the intensity of the CO(2-l) line to be anoma- 
lously low, i.e. below the 1:1 relation. These cases could be due 
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Figure 2. CO(l-O) and CO(2-l) IRAM 30 m spectra of the galaxies detected. The spectra have been binned to 30.6 km s — 1 and the scale is T m t, in mK. The 
CO(1-0) spectra are in the top half of eac h panel and the CO(2-l) spectra are in the bottom half. The arrow indicates the systemic velocity of stellar absorption 
lines, taken from our new SAURON data l lCappellari etal.l201lj) or from LEDA (for the 7 galaxies observed in CO but dropped from the ATLAS 3D sample). 
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Figure 2 - continued 
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to either subthermal excitation or pointing errors which cause the 
intensity of the 2-1 line to be underestimated. 

Molecular gas in nearby spirals is observed to have surface 
densities ranging from > 10 3 Mq pc -2 i n some galaxy nucl ei to 
a few tens of Mq pc -2 in the outer discs I IRegan et alj|2006l) . In- 
deed, sensitive observations of both HI and CO emission in spirals 



suggest that a typical cold gas disc in a spiral becomes dominated 
by molecular gas whe n the total gas su rface density is greater than 
about 14 Mq pc -2 jLerov et"aL 120081) . For comparison, Figure |4] 
shows the distribution of the average molecular surface densities 
among the CO detections of our sample. Since our data consist of 
a single pointing only, those surface densities are averages over the 
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30m beam (2 to 4 kpc in diameter), based on the intensity of the 1-0 
line. Mass densities are multiplied by a factor of 1.36 to account for 
helium. The molecular surface densities in the CO-rich ATLAS 3D 
galaxies are well within the range occupied by nearby spirals, and 
25% of our detections have average surface densities greater than 
50 Mq pc -2 . Galaxies showing signs of double-horned or asym- 



metric profiles (Table [5} have a surface density distribution which 
resembles that of the whole sample, with the caveat that such signs 
cannot be identified if the signal-to-noise ratio is low. The appar- 
ent dearth of double-horned or asymmetric profiles at small surface 
densities may only be an effect of the noise level. These considera- 
tions suggest that the CO discs in early-type galaxies probably have 
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gross characteristics which are not too different from the molecular 
discs in spirals. 

The galaxies with the highest average surface densities are 
NGC 1266 (230 M Q pc -2 ) and NGC4710 (210 M pc -2 ). Since 
the line ratios for NGC 1266 already indicate that the gas distribu- 
tion is more compact than the 1-0 beam, this galaxy's true molecular 



surface density is even higher than 230 Mq pc~ . In fact, high res- 
olution interferometric observations show the molecular gas in the 
galaxy is so co mpact that the surf ace densities are on the order of 
10 4 M pc" 2 dAlataloetal .11201 lh . 
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Figure 4. Average molecular surface densities for the detected galaxies. The 
open histogram shows all detected galaxies. Horizontal and vertical lines 
show galaxies with double-horned and asymmetric profiles, respectively. 



4.3 CO linewidths 

As mentioned above, the CO line shapes and widths carry infor- 
mation about the gas kinematics and radial extent, and appropri- 
ate analysis of the linewidth distribution can help extract this in- 
formation at least in a statistical sense. The 1300 km s" 1 band- 
width should be large enough to detect all reasonable CO lines, if 
the signal-to-noise ratio were not a limitation. In practice, of course, 
for a fixed H2 mass, spreading the emission over a larger velocity 
range will decrease the amplitude of the peak and make the line 
more difficult to detect. And since a broad CO line implies a mas- 
sive host galaxy with a deep potential well, a bias against the detec- 
tion of broad lines could potentially affect the statistical analysis of 
molecular mass distributions. 

In this context it is useful to consider the effects of our analy- 
sis procedures on the detectability of faint emission lines. A line of 
low amplitude may erroneously be included in the estimate of the 
spectral baseline level. Then, in searching for detections, we sum 
the baseline-subtracted spectrum over a 300 km s _1 velocity range, 
which (of course) will not be the entire line if the linewidth is larger 
than 300 km s _1 . Both of these procedures make the recovered line 
area systematically smaller than the "true" line area for weak lines. 
However, for lines of 300 km s _1 or narrower, the recovered area 
is a constant fraction of the "true" area regardless of the linewidth. 
Monte Carlo simulations of our analysis procedure confirm this ex- 
pected result. The simulations also confirm the strong bias against 
broader lines (at a fixed line area) for linewidths greater than 300 
km s _1 . For example, typical nondetections in our sample have 
thermal noise levels such that an integral over 300 km s" 1 has an 
uncertainty of 0.33 K km s _1 . A line of 3.3 K km s _1 area and 
600 km s _1 width (amplitude 5.5 mK), and a typical thermal noise 
level, actually has only a 40% chance of being recovered as a de- 
tection by our analysis procedures if it is mistakenly included in the 
baseline estimate, though it has a 97% chance if it isn't included in 
the baseline. Thus, for lines of 300 km s _1 or narrower there should 
be no bias in the detectability as a function of linewidth (for a fixed 
line area), though for lines broader than 300 km s^ 1 the effect can 
be severe at small amplitudes. 

Figure [5] shows the observed distribution of CO linewidths in 
the ATLAS 3D sample. As described in section 14.21 these are the 
FWHM of either a fitted Gaussian or a double-horned profile (in 
fact both methods give the same width for profiles which are ob- 
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vious double-horns). Three galaxies which each have two features 
detected at greater than 3a formal significance in a Gaussian fit are 
here interpreted as faint double-horned profiles; they are NGC4281 
and NGC 4643 in J=2-l (634 km s" 1 and 366 km s" 1 wide, respec- 
tivel y) and NGC4596 in J =l-0 (231 km s" 1 wid e). Six detections 
from lWelch & Sagel J2003I) and lSage et all < l2007t) do not have fitted 
linewidths quoted in their discovery papers, and in these cases the 
linewidths are estimated from plots of the spectra; four of the six 
are asymmetric enough that fits would be useless in any case. 

The observed CO linewidths can be compared to simple model 
distributions generated as follows. In "TF model 1," an estimated 
V c is generated for each ATLAS 3D galaxy from its Mk via th e 
SO Tully-Fisher relation of Iwilliams. Bureau & Cappellaril d20ld) . 
For each member of the sample we then generate 20 model CO 
linewidths as 2V C sin(i), where the inclinations are random and the 
probability of an inclination i is p(i) = sin(i). The model CO 
linewidth distribution is scaled down in amplitude by the factor of 
20 and by the global detection rate 0.22 to appear as the histogram 
in Figure [5] In "TF model 2," a similar procedure is followed ex- 
cept that V c is gener ated from the galaxy's dynamical mas s, again 
using the relations of lWilliams. Bureau & Ca ppellari (2010), and in 
this case only the CO detections are used. Additional details on the 
dynamical masses are given in Section |5T2l 

These models show the expected CO linewidth distributions 
under some non-trivial assumptions: (1) that the Tully-Fisher re- 
lations quoted above are appropriate; (2) that the CO is always in 
relaxed discs which reach to the flat, asymptotic value of the circu- 
lar velocity curve; and (3) in the case of model 1, that the CO-rich 
galaxies have the same Mk and circular velocity distribution as the 
entire sample. The third assumption is addressed in Section [5] be- 
low, where we show that the CO detection rate is not dependent 
on Mir. There is a modest bias in the CO detection rate with dy- 
namical mass, but model 2 already corrects for that bias by using 
only the CO detections. Some justification for the second assump- 
tion is suggested by the preponderance of double-horn profiles, but 
interferometric maps and independently-derived circular velocity 
profiles will be required in order to check for which galaxies the 
assumption is not satisfied. Four c ases are discussed in detail by 
lYoung. Bureau & Capp ellari (2008). 
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It is apparent that we have detected many more relatively nar- 
row lines, of FWHM 100-200 km s _1 , than would be expected in 
either of the simple models. Based on the above analysis of our 
detection procedures, we believe that the observed peak at around 
150 km s" 1 and the associated drop in the histogram between 150 
km s" 1 and 300 km s _1 are not imposed by our methodology. 
They could be related to a breakdown of the assumptions, which 
would mean that we don't have an accurate estimate of V c or the 
CO doesn't extend far enough in radius to trace V c in some cases. 
Using the profile full width at 20% of the peak, rather than at 50%, 
would eliminate some of the discrepancy but not all of it. For Gaus- 
sian profiles that change would increase the measured width by a 
factor of 1.5, but for double-horn profiles it should not affect the 
measured width. There are also several observational effects which 
would tend to produce an overabundance of narrow lines. Narrow 
lines will arise if the CO is extended and pointing errors prevent the 
detection of one horn of a double-horned profile (as in an extreme 
version of our CO 2-1 spectrum of NGC 3665). If the CO is in a 
ring whose radius is larger than the 30m beam, a narrow line may 
also result. Based on optical images, this effect may be a problem 
for a galaxy such as NGC 5379. Both of these discrepancies would 
also be rectified with interferometric maps. One impact of these ef- 
fects is that care is required when u sing the CO linewidths for a 
Tully-Fisher study jDavis etal.l201lh . 



5 MOLECULAR MASS CORRELATIONS 

One of the major unsolved questions concerning the molecular gas 
in early-type galaxies is its origin - whether it has been present in 
one form or another since the galaxies were assembled or whether 
it was more recently acquired from some external source. The 
ATLAS 3D data give a wealth of information on the structural, dy- 
namical, and stellar population properties of the galaxies, or in other 
words, their assembly and star formation histories. Examination of 
which types of galaxies are gas-rich and which are gas-poor may 
help to address the origins of the gas and, by implication, its role in 
the evolution of different kinds of early-type galaxies. Figures [6] 
through [14] compare the molecular gas contents to various other 
properties of the host galaxies. We discuss the implications of these 
results below. 



5.1 Stellar luminosity 

Figure|6]shows the distribution of molecular mass with if -band lu- 
minosity and a histogram of Mk for the CO-detected galaxies. CO 
emission is detected over nearly the entire range of luminosities in 
the sample, and there is no clear trend in M(H2) versus Mk - There 
is a small hint of a decrease in the detection rate at the highest lumi- 
nosities. Only one of the 13 galaxies with Mk < —25 is detected, 
whereas 3 would have been expected at the global average detection 
rate. However, the statistical significance of this difference is low. A 
stellar mass of 10 11 Mq corresponds to Mk = —24.4, at a typical 
stellar mass-to-light ratio of M*/L K = 0.82 dBell et al.ll2003h : the 
CO detection rate among galaxies brighter than this is 6/25 = 0.24 
±0.09, consistent with the global average. 

A Kolmogorov-Smirnov (KS) test on the Mk distributions of 
CO detections and nondetections indicates that they are consistent 
with each other (the probability that those two Mk distributions 
could have been drawn from the same parent population is 0.52). 
Thus we infer a constant CO detection rate at all luminosities from 



Mk = —21.5 to —26. Refuting that model will require observa- 
tions of many more galaxies brighter than Mk = —25, but such 
galaxies are rare in the local universe. 

The dependence of molecular mass on galaxy luminosity can 
also be investigated by comparing the cumulative H2 distribution 
functions for high, medium, and low luminosity galaxies (Fig- 
ure |7). Since the bulk of the CO observations produced non- 
detections ("censored" data) , the appropriate statistical tool is 
the Kaplan-Meier estimator for the cumulative distribution func- 
tion of a randomly censored sample. We calculate the Kaplan- 
Mei er estimator using the s o ftware packages ASURV version 
1.3 dlsobe & Feigelsorj 1 199(1 iFeigelson & Nelson] 1 19851) and R 
dR Development Core T eam 2009), with consistent answers from 
both packages. Here we divide the ATLAS 3D sample approxi- 
mately into thirds by luminosity, i.e. Mk ^ —23.4, —23.4 < 
M K < -22.5, and -22.5 < Mk- The CO detection rate in these 
bins is 19/92, 17/88, and 20/79, respectively. Figure|7]shows that the 
M(H2) distributions for these three luminosity bins are all consis- 
tent with each other, an impression that is confirmed by the Gehan, 
logrank, and Peto two-sample tests. Specifically, in the assump- 
tion that the H2 mass distributions for the bright and faint galax- 
ies are drawn from the same underlying population, the probability 
of measuring a difference as large as the one observed is estimated 
at 0.57 to 0.6 0. For comparison, the "dichotomy" stellar mass of 
3 x 10 10 M dKauffmann et alj2003l) occurs at M K ~ -23.1, for 
M- t / Lk = 0.82 (Bell et al. 2003); this stellar mass is in the middle 
luminosity bin. Our analysis indicates no measurable difference in 
the H2 content of galaxies as a function of their stellar luminosity. 

The use of the Kaplan-Meier methodology merits some dis- 
cussion, as it assum es the censoring pattern is random (e.g. 
IWall & Jenk ins 2003). In this case the censoring in brightness tem- 
perature is certainly not random, as observations were made to a 
fixed brightness temperature noise level. However the censoring in 
H2 mass and M(H2)/ Lk should be approximate l y random, for the 
reason s dis cussed by llsobe & Feigelsorj dl992l) . lO'Sullivan et alj 
( l200lh . and lWall & Jen kins (2003). Multiplication by the square of 
distance and division by the optical luminosity effectively serve to 
randomize the censoring pattern in those derived quantities, since 
our sample spans a range of distances and there is no hint of a 
dependence of CO emission on optical luminosity in our sample 
(Figure [6}. Thus we believe the assumption of random censoring 
is justified. In addition, our objectives are to assess the similarity 
of the cumulative distribution functions of different subsets of our 
data and there is no reason to believe that the censoring should dif- 
fer from one subset to another. 

As mentioned above, the practical difficulty of detecting a 
broad and weak spectral line means that a given mass of molec- 
ular gas will be easier to detect in a face-on galaxy or in an edge- 
on galaxy with a small circular velocity than in an edge-on galaxy 
with a circular velocity of 200 km s _1 or higher. This effect could 
cause the H2 mass distribution to be underestimated in the highest 
luminosity bin. Careful simulations would be required to gauge the 
magnitude of the effect, but the sense of the correction would be 
to increase the incidence of molecular gas in the highest luminosity 
galaxies. 

The lack of a correlation between H2 mass and stellar lumi- 
nosity has been interpreted in the past to mean t hat the molecular 
gas is unrelated to interna l stellar mass loss (e.g. Knapp & Rupen 
1996; IWiklindetal.i ri995). The assumption is that internal mass 
loss would produce a linear relationship between stellar luminosity 
and H2 mass. However, if the mass loss material is shock-heated 
to 10 6 K or higher by the relative stellar velocities or is other- 



© 0000 RAS, MNRAS 000, 000-000 



Molecular gas content of early-type galaxies 



13 



9.5 
9.0 
8.5 
8.0 
7.5 
7.0 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



5 •••••• . . 

ft*:£»V ********* 

% * ** + ™ + * * ♦ » 

.•;•.>:•«*»•.!:••.• * 

+ ++■*■ + 

+ + + 

+ 

" I I I I I" 



-21 



■22 -23 -24 -25 



■26 



60 

50 

40 

30 

20 

10 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 III 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-21 



' '''I I I I I ill ill I |>1 



14 

12 

10 

8 

6 

4 

2 




-22 -23 -24 -25 



-26 



M K M K 

Figure 6. Molecular masses and Mj(. Small crosses are ATLAS 3D sample galaxies that are not detected in CO emission (3<r upper limits), and large circles 
are detections. The histograms in the right-hand panel show both the properties of the entire ATLAS 3D sample (dotted line; left-side scale) and those of the 
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Figure 7. The Kaplan-Meier estimators for the cumulative M(H2) distri- 
bution functions (Sk) for high-, medium-, and low-luminosity ATLAS 3D 
galaxies. These estimators include the effects of censored data (upper limits) 
in constructing the cumulative distribution function. The shaded regions in- 
dicate 1<t uncertainties on the estimated cumulative distribution functions. 
The distribution functions are similar for all three luminosity ranges, sug- 
gesting no dependence of molecular mass on host luminosity. 



wise heated by ambient hot gas, it would be difficult to predict 
how much of it might be able to cool and re-form molecules. Envi- 
ronmental and feedback effects would also complicate the relation- 
ship between the stellar luminosity and the amount of molecular 
gas present in a galaxy. Thus, we argue that one cannot infer either 
internal or external origin for the gas based on the M(H2)-Mif re- 
lationship alone (or the lack of one). The value of the ATLAS 3D 
project is, of course, that other structural and kinematic parameters 
and stellar population information are available for the sample as 
well. 

5.2 Dynamical mass 

In spite of the lack of a relationship between CO content and Mk, 
there is some evidence for a modest dependence of the CO content 
on the dynamical mass and a closely related quantity, the global 



stellar velocity width a e . Figure [8]shows that the CO detection rate 
for galaxies with a e ^ 100 km s _1 is relatively high, 0.38 ± 0.05, 
whereas for a e > 100 km s _1 it is 0.14 ± 0.03. The median a e 
for CO detections is 99 km s _1 , whereas for nondetections it is 
137 km s" 1 . The effect in dynamical mass is more subtle, so that 
the CO detection rate above a dynamical mass of 3 x 10 10 M is 
0.19 ± 0.03 and below is 0.25 ± 0.04. A KS test on the dynamical 
masses of CO detections and nondetections gives a probability of 
6.5% that they could have been drawn from the same parent pop- 
ulation. However, the Kaplan-Meier estimators show no significant 
differences in H2 masses of high and low mass galaxies. 

For both o c and dynamical mass the sense of the difference is 
that lower mass galaxies are more likely to be detected. It is curious 
that the magnitude of the effect is much stronger in <r e than in dy- 
namical mass. The details of the scaling relations between those two 
quantities are beyond the scop e of this paper, but it is helpful to be 
precise about their definitions. Cappe llari et al.l (|2006) explain that 
<T e is the second moment of the galaxy's luminosity-weighted line 
of sight velocity distribution (LOS VD). Because it is constructed by 
fitting a Gaussian to a stack of all spectra within the effective radius, 
it includes contributions from both the mean stellar rotation and the 
velocity dispersion approximately as a e ~ {V^ ot + a 2 ) 1 / 2 , where 
Vrot is the projected mean st ellar rotation velocity an d a is the lo- 
cal stellar velocity dispersion. ICappellari etail ( 2006 ) discusses the 
use of o e in the virial mass as M V i r ~ 5a 2 i? e /G- The dynamical 
mass used here is, however, derived fr om stellar dynamical models 
of the SAURON stellar kinematics dCappellari et al, 1 12010 ), using 
the Jeans anisotropic modelling technique of ICappellaril (120081) . It 
can be understood as A/jam ~ 2 x M1/2, where M1/2 is the total 
dynamical mass within a sphere containing half of the galaxy light. 
In addition, as M1/2 is generally dominated by the stellar mass, 
Mjam approximates the stellar mass M ata r but should not be con- 
flated with the total halo mass out to the virial radius. The quantities 
Mjam and a e for the ATLAS 3D sample will be published in Cap- 
pellari et al., in preparation. 

Several different effects could be responsible for the depen- 
dence of the CO detection rate on a e and dynamical mass. We note 
first that molecular gas likely does induce star formation, and the 
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younger stellar populations should be both brighter and dynami- 
cally colder than the rest of the stars. Star formation should not af- 
fect the dynamical mass or cr e , though, as <r e already includes both 
the rotation and the dispersion components of the velocity width. In- 
stead, the dependence of the CO detection rate on dynamical mass 
and a e could be driven by either or both of (1) an observational 
bias against the detection of broad CO lines, as described above, 
and (2) a downsizing or rejuvenation effect in which the more mas- 
sive galaxies are less likely to contain molecular gas. The obser- 
vational bias would manifest itself as a deficit of CO detections in 
edge-on high mass galaxies compared to face-on high mass galax- 
ies. A downsizing effect would be consistent with much other recent 
work on star formation histories as a function of galaxy mass (e.g. 
Rogers et al.l2010l : Izhu, Blanton & Moustakashoiol : iThomas et all 
20101). 

Of course, regardless of whether the dependence on dynamical 
mass is an observational bias or a real downsizing effect (or both), 
something of a 'conspiracy' between the mass and the mass-to-light 
ratio is required in order to reproduce the result that the CO detec- 
tion rate is completely independent of Mk- This issue can also be 
investigated in greater detail with estimates of the stellar population 
mass-to-light ratios of the ATLAS 3D galaxies. 

5.3 Morphological type 

There is a significant difference in the CO detection rate between 
the galaxies traditionally classified as elliptical and lenticular. The 
ATLAS 3D sample contains 68 galaxies whose morphological type 
in LEDA is T < -3.5 (ellipticals). Only 3 of the 56 total CO de- 
tections are in this group, for a detection rate of 0.05 ± 0.03. These 
detections are NGC2768, NGC 5173 and N GC 4283, and the lat- 
ter is of marginal quality dSage et al.ll2007h . However, we caution 
that these morphological types are generally based only on photo- 
graphic imaging and have the limitations (or perhaps "features") 
that they carry little kinematic information and that the presence of 
a dust lane in the optical photographs could cause a galaxy to be 
classified as lenticular rather than elliptical. 

Interferometric imaging shows that the molecula r gas in early- 
type galaxies is invariably associated wi th dust lanes (Youngl l2002l . 
l2005tlYoung. Bureau & Cappellaril2008h . so that if the morpholog- 
ical classification is indeed biased by the presence of a dust lane 
then it is nearly a foregone conclusion to find a higher CO detec- 
tion rate in lenticulars than in ellipticals. Indeed, Paper II tabulates 
the evidence for dust discs, dusty filaments, and blue regions in op- 
tical images of the ATLAS 3D galaxies. Of the galaxies with one 
or more of such features, 3/49 (6%) are classified as elliptical and 
46/49 as lenticular; but of the galaxies with no dust or blue fea- 
tures, 64/209 (31%) are elliptical and 145/209 are lenticular. Of 
course these statistics do not prove a dust- or star formation-bias 
in morphological classification, but they are consistent with such an 
interpretation. 

From a different point of view, the SAURON observations of 
kinematics in early-type galaxies have shown that there is often lit- 
tle if any structural di fference between lenticulars and ellipticals 
temsellem et al.l [2004). so a purely isophotal classification con- 
tributes little to our understanding of their formation histories. 



5.4 Specific angular momentum 

lEmsellem et alj J2007l) advocate that a more fundamental way to 
classify early-type galaxies is with the Xr parameter, a simplified 



and dimensionless version of the luminosity weighted specific an- 
gular momentum. It encapsulates some of the information on the 
degree of ordered vs. random motions (i.e. rotational vs. pressure 
support) that one gets from the classic V ro t/o- ratio. In addition, as 
some systems show rotation only within a few hundred pc of the 
nucleus, the Xr parameter also reveals something about the spa- 
tial extent of the rotation through the radial contribution to the spe- 
cific angular momentum. Slow rotators (including galaxies with ro- 
tation only at small radii) have small values of Xr. In an update 
of this earlier work, Paper III demarcates fast and slow rotators 
at Xr = O.Sly'e, where e is the photometric ellipticity and both 
quantities are measured as luminosity-weighted averages over the 
effective radius. Paper III also discusses the relationship between 
the traditional E/S0 and the slow/fast rotator classifications. There 
is some degree of alignment between the systems, as the slow ro- 
tators are mostly classified as ellipticals. However, the converse is 
not true, as 66% of the ellipticals in ATLAS 3D are fast rotators. The 
consistency between the systems is only strong at the extremes of 
Xr and e. 

Figure|9]shows normalized H2 masses and CO detection rates 
as a function of the ratio Xr/ '(0.31-^/e), and Table Q] shows that 
the CO detection rate among slow rotators (0.06±0.04) is signif- 
icantly lower than that among fast rotators (0.24±0.03). The CO 
detections in slow rotators are NGC 1222 and NGC 4476. Of these 
two, NGC 1222 is kinematically disturbed due to a recent interac- 
tion (Paper II). NGC 4476 also presents a complex structure as its 
stellar velocities drop towards the outer regions of the SAURON 
field and show signs of reversing at yet larger radii (A. Crocker, pri- 
vate communication). In short, there are few detections among the 
slow rotator class. 

Both a galaxy's value of Xr and its projected circular velocity 
V c sini are affected by inclination, of course. However, if a galaxy 
has small Xr due to being face-on, and if it has molecular gas, it 
will have a relatively narrow CO line and will be easy to detect. 
Thus the deficit of CO detections among galaxies of small Xr is 
not due to inclination effects because they have the wrong sense to 
reproduce such a deficit. 

The deficit of CO detections among galaxies of small Xr is 
also not driven by the trend with <r e or the dynamical mass. Pa- 
per III shows that the slow rotators do tend to be more massive 
than fast rotators, and they have higher <r e values, but the overlap 
between fast and slow rotators is sufficient that this effect can be 
controlled. Of the 36 slow rotators in the ATLAS 3D sample, we be- 
gin with the subset of 31 which have o e < 260 km s _1 or the 34 
which have 9.5 < logMjAM/M0 < 11.65. These subsets have 
CO detection rates of 0.06 ± 0.04, consistent with Table[T] We con- 
struct a matched sample of fast rotators by selecting (at random) 
one fast rotator which has <j e within ±12 km s _1 or Mjam within 
0. 1 dex of each slow rotator. By design, then, the matched sample of 
fast rotators has a distribution of Mjam or a e which is statistically 
indistinguishable from that of the slow rotators. Repeating the pro- 
cess 10 4 times, we find that the CO detection rate in our matched 
samples of fast rotators is 0.24 ± 0.07 when the matching is done 
by o~ e or 0.36 ± 0.07 when matched by A/jam- Here quoted un- 
certainties are dispersions in the set of 10 4 trials. In short, the CO 
detection rate is still significantly higher among fast rotators than 
among slow rotators, even after controlling for the a e and Mjam 
distributions. 
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Figure 8. Molecular mass, dynamical mass, and velocity width a e . Symbols and lines are as in Figure[6] Dynamical masses are obtained from Jeans anisotropic 
models (Cappellari et al. in prep) as described in Cappellari 1 2008). 
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5.5 Internal kinematic structure 

Paper II documents photometric evidence for internal substruc- 
tures such as bars, rings, shells, tidal tails, and anemic spiral struc- 
tures in the ATLAS 3D sample. We find no evidence for a differ- 
ence in CO detection rate between galaxies showing these types of 
features and galaxies without them (Table [Q. Similarly, we mea- 
sure the misalignment angle between the stellar photometric and 
kinematic axes, as these misalignments can be indications of non- 
axisymmetric components like bars. A KS test reveals no significant 
difference in the misalignment angles of CO detections and nonde- 
tections. 

Paper II also uses the features in the stellar velocity maps to 
classify the internal kinematic structure of galaxies. Galaxies are 
primarily classified as "regular rotators," or RR, when the velocity 
field resembles that of a thin, inclined disc such as a spiral galaxy. 
Others are "non-regular," or NRR. The RR class has a nearly one- 
to-one overlap with the fast rotators as defined by ellipticity and 
Xr, so these are assumed to be the same family of galaxies. Table 
Q]shows that the CO detection rate is significantly higher in the RR 
class than in the NRR class, and the 'matched sample' technique 
described in the previous section again confirms this result after 
controlling for different <r e and dynamical mass distributions of the 
RR and NRR classes. 

In terms of kinematic structure, early-type galaxies are quite 
heterogeneous (Paper II); they include galaxies with no observ- 
able rotation at all (non-rotators like M87), galaxies with kinemat- 
ically distinct or even counterrotating cores, galaxies comprised 
of two counterrotating stellar discs, such as IC0719, and (infre- 
quently) face-on disc-like regular rotators. Table [T] shows no par- 
ticularly strong evidence that the CO detection rate differs among 
these subtypes. The CO detections in NRR galaxies are NGC 1222, 
NGC 3073, IC 0719, and NGC 7465, of which NGC 1222 is also a 
slow rotator and is discussed in the previous section. Curiously, the 
"two a peak" galaxies made of two counterrotating discs are ex- 
ceptions to the general rule that the CO detection rate is higher at 
low cr e . This type has the lowest dispersions o e of any of the kine- 
matic classes. Their relatively low CO detection rate (compared to 
the low <r e galaxies in general) should not be an effect of the obser- 
vational bias against broad lines because, given their small masses, 
they should not have broad CO lines. 

Papers II and III clearly show that early-type galaxies have 
a variety of different assembly histories. Some general patterns of 
assembly are conducive to retaining a high specific angular momen- 
tum within an effective radius whereas other patterns are not. The 
patterns which develop low specific angular momentum can include 
misaligned major mergers, or mergers in which the orbital angular 
momentum cancels the spin angular momenta of the progenitors 
(Bois et al. in prep), or a large number of randomly oriented minor 
mergers. Our analysis suggests that galaxies with these types of as- 
sembly histories are least likely to be detected in CO. These may 
be the cases in which the misaligned angular momenta of gaseous 
discs coming into the merger caused the gas to drop to the center of 
the galaxy where it was consumed or destroyed. Alternatively, they 
may have lost their cold gas prior to their assembly or they may 
be less likely to acquire cold gas after assembly. Future work on 
gas kinematics and star formation histories can help to distinguish 
between these scenarios. 



Table 1. CO detection rates for kinematic types 



Type CO no CO Total CO det. rate 



Notes: Kinematic classifications are taken from Papers II and III. This table 
and the next exclude NGC 4486A, for which we have no CO data. 

5.6 Environment: Virgo Cluster 

The ATLAS 3D sample contains 58 galaxies that are within 3.5 Mpc 
of M87 and are considered members of the Virgo Cluster. Many 
of these galaxies (including M87) h ave distance measurements 
from surface brightness fluct uations dTonrv et al .11200 it iMei et all 
120071 ; ICappellari et al.11201 lab . When a surface brightness distance 
is available, it is used in the computation of a galaxy's distance to 
M87. Where independent distance measurements are not available 
for known Virgo Cluster members, a distance of 16.5 Mpc is used as 
that i s the adopted mean distance of subclusters A and B dMei et al] 
120071) . For comparison, those authors also quote the surface bright- 
ness fluctuation distance to M87 as 17.2 Mpc, the typical distance 
uncertainties for their data are 0.7 Mpc, and the line-of-sight depth 
of the cluster is 0.6 Mpc ±0.1 Mpc (lcr). We note also that there 
is no statistically significant difference in the Mk, dynamical mass, 
or a e distributions of cluster and non-cluster galaxies in ATLAS 3D , 
so there is no detectable mass segregation in this sample. 

Twelve of the 57 Virgo Cluster members with CO data are 
detected, for a cluster detection rate of 0.21 ±0.05 or less than lcr 
lower than the detection rate for the sample as a whole. However, 
since the Virgo Cluster is relatively nearby and we have observed 
our sample to a fixed noise level, it is more appropriate to compare 
the detection rate in Virgo to that of galaxies at a similar distance. 
The entire ATLAS 3D sample contains 89 galaxies at D ^ 20 Mpc 
and 124 galaxies at D < 24 Mpc; thus, within 24 Mpc, half of 
the sample galaxies are Virgo Cluster members and half are not. A 
restricted comparison of galaxies within 24 Mpc mitigates the ef- 
fect of distance, since the typical H2 mass limit at 24 Mpc is only a 
factor of 2 higher than at 17 Mpc. Table[2]compares the CO detec- 
tion statistics inside and outside of the Virgo Cluster for all galax- 
ies within 24 Mpc, and also for the subset of fast rotators within 
24 Mpc. The detection rate among cluster members (0.21 ±0.05) 
is nominally lower than the detection rate of non-members (0.29 
±0.06) and the detection rate of the entire sample (0.22 ±0.03). 
However, the difference is at only the la level in the combined 
uncertainties so more detailed investigations of the CO luminosity 
functions are appropriate. 

The top two rows of Figure [Tolpresent the normalized masses 
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0.20±0.03 




Kinematic subtypes 




Non-rotator 


1 


6 


7 


0.14±0.13 


KDC/CRC 


1 


18 


19 


0.05±0.05 


2 a peaks 


2 


9 


11 


0.18±0.11 


NRR/no feature 


1 


11 


12 


0.08±0.08 




RR subtypes 






RR/2max 


11 


25 


36 


0.30±0.08 


RR/other 


40 


137 


177 


0.23±0.03 
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Table 2. CO detection rates for galaxies within 24 Mpc 





CO 


no CO 


Total 


CO det. rate 


Virgo Cluster 


12 


45 


57 


ft 9 1 4-fl OS 


non-cluster 


19 


47 


66 


0.29±0.06 


Total 


31 


92 


123 


0.25±0.04 




Fast rotators only 




Virgo Cluster 


11 


37 


48 


0.23±0.06 


non-cluster 


19 


42 


61 


0.31 ±0.06 


Total 


30 


79 


109 


0.27±0.04 



M(H 2 )/L/c versus the local bright galaxy density, both for the en- 
tire ATLAS 3D sample and for the subset of galaxies closer than 
24 Mpc. Normalized H2 masses are also plotted as a function of 
the deprojected distance to M87. Paper VII explains the calculation 
of the local density pro as the number density of galaxies brighter 
than Ma- = — 21.5, averaged within a sphere which contains 10 such 
neighbors and quoted in units of Mpc -3 . In the ATLAS 3D parent 
sample, which includes all morphological types, all galaxies outside 
the Virgo Cluster have densities log pio < —0.4 and all but two of 
the Virgo Cluster members have log pio > —0.4. A KS test on the 
distributions of local densities for CO detections and nondetections 
(both within 24 Mpc; middle row, right side panel of Figure |10t 
gives a 31% probability that they could be drawn from the same 
parent population. 

The preceding evidence shows that the CO detection rate is 
nearly the same for early-type galaxies in the Virgo Cluster and 
outside the cluster. Clearly the galaxies in the cluster are not en- 
tirely devoid of molecular gas. However, it is still possible that they 
may have suffered some modest stripping, which would show up 
as systematically smaller H2 masses or M(H2)/La ratios for clus- 
ter members than for non-members. In order to test this hypothe- 
sis, Figure [TT] presents the Kaplan-Meier estimators for the cumu- 
lative H2 and M(H2)/La- distribution functions for galaxies inside 
the Virgo Cluster, galaxies outside the cluster but closer than 24 
Mpc, and the entire ATLAS 3D sample. In these analyses the null 
hypothesis is that the cluster and non-cluster galaxies have the same 
underlying M(H2) and M(H2)/La- distributions, and a test statistic 
as large as the observed one would be expected 25% to 60% of the 
time for M(H 2 ) and 46% to 57% of the time for M(H 2 )/L A '. These 
probabilities are not so small that the null hypothesis can be con- 
fidently excluded. In other words, the available data do not show a 
striking or reliable difference between the M(H2) and M(H2)/La 
distribution functions for cluster and non-cluster galaxies. 

While there may be a slight decrease in the CO detection rate 
in the Virgo Cluster compared to field early-type galaxies within 
24 Mpc, there is as yet no good evidence that the molecular gas 
contents are systematically smaller inside the cluster. This result 
stands in opposition to the situation for atomic gas; HI detection 
rates for early-type galaxies in the Virgo Cluster are acutely low 
compared to detection rates in the field ( di Serego Alighieri et alj 
120071 : iGrossi et al.ll2009l; lOosterloo et ai1l2010l) . Additional discus- 
sion of galaxy evolution in the cluster can be found in Section[6] 

Figure [T2] shows the systemic velocities of CO detections and 
non-detections as a function of projected distance to M87, as well 
as the velocity distributions of ATLAS 3D Virgo Cluster members 
and CO detections. The Virgo Cluster members have a centrally 



peaked distribution with a velocity dispersion of 525 km s~ , in 
good agreement with the value of 573 km s" 1 measured for el- 
liptical and lenticular galaxies by Bin ggelli. Tammann. & Sandagej 
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Figure 11. The Kaplan-Meier estimators for the cumulative M(H2) and 
M(H 2 )/La distribution functions (S K ) for ATLAS 3D Virgo Cluster mem- 
bers, galaxies outside of the Virgo Cluster (but at distances D ^ 24 Mpc 
from Earth), and for all ATLAS 3D members. The shaded regions indicate 
lcr uncertainties on the estimated distribution functions. 



Jl987h . The 12 Virgo Cluster CO detections have a dispersion of 
414 km s _1 and are also peaked about the cluster systemic veloc- 
ity. This behavior is in marked contrast to the broad, flat velocity 
distribution of the spirals in the cluster, which are characterized 
by a dispersion of 888 km s _1 (Binggelli, Tam mann. & Sa ndage 
1 1987). A more compl ete characterization of the spirals is provided 
by Rubin et aT] J 1999h . who find that the disturbed spirals are ap- 
proaching a relaxed velocity distribution whereas the regular, undis- 
turbed spirals have a broad, flat, non-virialized velocity distribution. 

The dynamical status of the CO-detected Virgo Cluster early- 
type galaxies can be further assessed by comparing their veloc- 
ity distribution to those of the other early-type galaxies and the 
late-type galaxies in the cluster. A comparison sample of late- 
type galaxies in the Virgo Cluster is provided by the Sa-Sm, Im, 
BCD, Sp, and dS galaxies in the Vir go Cluster Catalog (VCC; 
Binggelli, Sandage. & Tammann! 1 1985h . For this purpose we se- 
lect from the VCC only the galaxies identified as cluster members 
and having measured systemic velocities. The corresponding early- 
types are the VCC E, SO, dE, and dSO. Velocity distributions for 
these samples are shown in Figure [13] The results of KS tests on 
the velocity distributions of these various samples are shown in 
Table [3] Notably, the velocity distribution of the ATLAS 3D Virgo 
Cluster CO detections is consistent with those of the ATLAS 3D 
Virgo Cluster CO non-detections and the VCC early-types. How- 
ever, our ATLAS 3D cluster CO detections are not consistent with 
the VCC late-type galaxies; that null hypothesis can be rejected 
at the 95% confidence level. Thus, while the number of CO de- 
tections among the ATLAS 3D Virgo Cluster galaxies is still rela- 
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Figure 10. Molecular gas properties inside and outside the Virgo Cluster. The quantity pm is a local galaxy density presented in Paper VII, and the dashed 
line at logpio = —0.4 separates Virgo Cluster members and nonmembers (Paper I). Other symbols and histograms are as in Figure [6] The middle row is 
restricted to galaxies with distances less than 24 Mpc (mitigating the effects of distance on CO detection). The third row presents molecular masses against the 
three-dimensional (deprojected) distance to M87. 
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Figure 12. Systemic velocities of CO-detected and non-detected galaxies near the Virgo Cluster. Symbols in the right-hand plot are as in Figure |6] The right 
panel shows systemic velocities of all ATLAS 3D Virgo cluster members in the open histogram (heavy line) and systemic velocities of the CO detections in the 
hashed histogram. Unlike the previous histograms, in this case the vertical axes for detections and non-detections are not scaled relative to each other. 
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Figure 13. Velocity distributions for ATLAS 3D Virgo Cluster members, 
VCC early-type galaxies, and VCC late-type galaxies. The ATLAS 3D CO 
detections are indicated in the hashed histogram. ATLAS 3D galaxies are 
binned to 200 km s~ 1 ; VCC galaxies have smaller bins whose sizes are ad- 
justed so that their histograms occupy the same area as that of the ATLAS 3D 
Virgo Cluster members. 



Table 3. Kolmogorov-Smirnov tests 



Sample 


A3D Virgo Cluster 


A3D Virgo Cluster 




(all) 


(CO detections) 


VCC late-types 


0.021 


0.052 


VCC early-types 


0.863 


0.262 


A3D CO nondet. 


0.999 


0.228 



Notes — 



Matrix of probability values for the KS test statistic on the systemic veloc- 
ity distributions of ATLAS 3D galaxies and members of the Virgo Cluster 
Catalog. 

place them within a true distance of 1.3 Mpc from M87. Thus, they 
must be physically located within the hot intracluster medium. In 
combination with the centrally peaked velocity distributions, these 
data suggest that the CO detections in the Virgo Cluster are indeed a 
dynamically relaxed population that does not avoid the center of the 
cluster. This behaviour is markedly different to that of atomic gas 
in the spirals of the Virgo Cluster, which show strong HI deficien- 
cies n ear the center of the Cluster dCavatte et al.ll990tlChung et alj 
l2009t) . 



tively small, the indications are that the dynamical properties of the 
CO-rich early-type galaxies are more consistent with the virialized 
early-type galaxies than with the spirals and irregulars. It thus ap- 
pears unlikely that our CO detections are in new or recently added 
cluster members. 

Figure [14] shows the spatial distribution of ATLAS 3D targets 
and CO detectio ns in the Virgo cluste r, compared to the location 
of the hot gas jBohringer et al.|[l994h . Several of the CO detec- 
tions are seen in projection in the central regions occupied by hot 
gas. Notable among these are NGC4429, NGC4476, NGC4435, 
NGC 4477, and NGC 4459, all of which are secure CO detections 
and are within 2° (0.6 Mpc) of M87. NGC 4435 and NGC 4477 
are members of the Markarian chain to the northwest of M87. 
NGC 4435, NGC 4476, and NGC 4459 have distance me asurements 
from surface brightness fluctuations (IMei et al.l 12007). and these 



5.7 Groups and low density regions 

Figure [Tol shows possible evidence for an enhanced CO detection 
rate at the outskirts of the Virgo Cluster, « 15-20° (5-8 Mpc) from 
M87. This result is driven by the galaxies near NGC 4278 (north 
of M87), the Leo II or NGC 3607 group (northwest of M87), and 
galaxies near NGC 4684 (south of M87). For the galaxies 5 to 8 
Mpc from M87 the CO detection rate is 8/17 = 0.47 ± 0.12, not 
quite 2cr above the global detection rate. 

Paper VII also presents measurements of a local density E3 for 
the ATLAS 3D members; this estimate gives the projected density of 
galaxies within a cylinder of height 600 km s _1 and radius such that 
3 neighbors are inside. As discussed in Paper VII, E3 is large both 
for members of the Virgo Cluster and for galaxies in small, tight 
groups. A KS test on the E3 distributions of CO detections and 
nondetections (at distances less than 24 Mpc) gives a probability of 
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RA (J2000) 

Figure 14. CO detections and non-detections towards the centre of the 
Virgo Cluster. The gre yscale and contours are the ROSAT X-ray image of 
iBbhringer et al.1 l ll994l) . The red crosses are the ATLAS 3D non-detections 
and blue filled circles are the CO detections. The CO detection nearest the 
X-ray peak is NGC 4476, whose regular and apparently undisturbed CO disc 
was mapped bv lYoundEioi . 



7% that the two distributions could have been drawn from the same 
parent - not a compelling difference, as we found for the cluster 
vs. non-cluster comparisons in the previous section. Future studies 
investigating the effect of the group environment on CO content or 
pre-processing of molecular gas in galaxies infalling into the cluster 
will require larger samples. 

A strong trend in CO detection rate with local density can only 
be found when considering the lowest density environments and 
lowest mass galaxies in this sample. The CO detection rate rises 
to 0.55 ± 0.09 for galaxies of a e < 100 km s _1 and density 
log(/9io/Mpc -3 ) < —1.7, or 0.41 ±0.11 for galaxies of dynam- 
ical mass < 2 x 10 10 Mq in the same density regime. However, 
as Figure [TO] suggests, the lowest density environments primarily 
occur at large distances in our sample and they drop out of the sub- 
samples limited to 24 Mpc. In addition, the trend with local den- 
sity is not significant for galaxies of a e > 100 km s -1 or mass 
> 2 x 10 10 Mq. Thus, this increase in the CO detection rate of 
low mass, relatively isolated early-type galaxies did not manifest 
itself in the previous discussions of environmental effects. 

The half-dozen most CO-rich galaxies of the sample do tend 
to be in relatively low density environments, which is undoubt- 
edly an important clue to their formation histories. This CO-rich 
subsample includes NGC 1266, NGC 2764, NGC 1222, NGC 3665, 
UGC 09519, and NGC 6014. These galaxies are all more distant 
than 20 Mpc, in regions where the local luminosity density is a 
factor of 10 to 200 lower than that in the Virgo cluster, and the 
distances to their nearest Mk < —21.5 neighbors are ^ 1 Mpc. 
If the accretion of cold gas is most efficient in these kinds of poor 
environments, such accretion might provide a natural explanation 
for the most CO-rich ATLAS 3D galaxies. Signatures of the accre- 
tion could be sought in the CO kinematics (through interferometric 
imaging) or in the gas-phase metallicity. 



5.8 Young stellar populations 

As one might expect, the molecular gas content of early-type galax- 
ies is tightly correlated with dust and young stellar populations, 
which are indicated for the ATLAS 3D sample in Paper II. The CO 
detection rate is 37/48 (0.77 ± 0.06) in galaxies with dust discs, 
dusty filaments, or blue regions, but only 17/208 (0.08 ± 0.02) in 
galaxies with none of those features. 

Figure \T5\ shows that the galaxies that are rich in molecular 
gas also have evidence of young stars, as traced by the H/3 absorp- 
tion line strength index in the classic Lick/IDS system from our 
SAURON spectroscopy (McDermid et al. in prep). Here the H/3 ab- 
sorption line strength is measured in an aperture of radius R e /2, 
which has a median value of 9.3" for the sample so is similar to the 
spatial region in which we have searched for CO J=l-0 emission. 
However, the trend is not qualitatively different when apertures of 
radius R e or R e /8 are used. Naturally, as H/3 is an equivalent width 
measure, both it and M(H2)/Lk are roughly normalized to the total 
stellar luminosity. 

There is a clear trend for galaxies with strong H/3 absorption 
to exhibit higher M(H2)/La _ ratios, indicating that higher gas frac- 
tions are found in galaxies with larger proportions of young stars. 
Intuitively, the trend makes sense, since the H2 reservoir is thought 
to fuel any ongoing star-formation. However, the relation is not es- 
pecially tight, with more than an order of magnitude variation in 
M(H2)/Ljf at a given line strength. This spread probably reflects the 
differences between these two tracers. The mass of H2 should trace 
the instantaneous star-formation activity, with some caveats noted 
below; H/3, in combination with metal indices, traces a weighted 
average age which is strongly biased towards recent star formation 
activity. Differences in recent star formation histories (especially 
bursty or merger-driven activity) could produce some of the scatter. 
Some scatter is also probably due to the fact that H/3 is degenerate 
for a broad range of possible star formation histories and does not 
exclusively trace recent (< 1 Gyr) events. 

Another cause of scatter in the H/3-H2 trend could be differing 
star formation efficiencies, driven by several mechanisms. For ex- 
ample, the shape and depth of the pote ntial can affect the stability 
of the cold gas disc embedded within it dMartig et al .1 2009). It may 
also be that the molecular gas reservoir is prevented from forming 
stars in some galaxies due to some form o f heating, e.g. from A GN 
activity or hot phases of stellar evolution. ICrocker et al.l d20ld) ex- 
plore the latter in some detail, proposing connections to the optical 
emission from the ionized gas. We postpone such analysis using the 
full ATLAS 3D data complement for future papers in this series. 

There are a number of CO upper limits found in galaxies with 
strong H/3 absorption (> 2.3 A). At face value, these facts could 
indicate the presence of a significant young stellar population with- 
out any associated molecular gas, perhaps immediately following 
exhaustion of the gas after a significant star-formation event. How- 
ever, those H2 upper limits are quite high because the galaxies are 
distant. These non-detections may be comparable in their H2 prop- 
erties to the detected galaxies, simply falling just below our sen- 
sitivity level. At the other extreme, while there are many objects 
with intermediate values of H/3 that are both CO-detected and un- 
detected, there are no CO detections in objects with H/3 1.4 
A. The galaxies showing no evidence for young stellar populations 
also show no evidence of molecular gas within our detection thresh- 
olds. 



© 0000 RAS, MNRAS 000, 000-000 



Molecular gas content of early-type galaxies 21 



10"' r 



10"" r 



10" 4 r 



iO _:> La 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



i i i 



12 3 4 5 

H/S abs. (A) R e /2 

Figure 15. Molecular mass and H/3 absorption line strength (McDermid et 
al. in prep). Symbols are as in Figure|6] 



6 DISCUSSION 
6.1 Rejuvenation 

iThomas et ail (2010) have used SDSS spectroscopy to measure 
luminosity-weighted simple stellar population (SSP) ages for a 
sample of low redshift, morphologically selected early-type galax- 
ies. They interpret a population of younger galaxies as examples of 
a rejuvenation phenomenon driven by recent (< 2 Gyr) star forma- 
tion. They also argue that this rejuvenation phenomenon is mass- 
and density-dependent, being more common (or perhaps stronger?) 
in low mass galaxies in low density environments. 

This picture is at least qualitatively consistent with what we 
find for the molecular gas content of early-type galaxies, with a few 
important caveats. For example, the results presented above suggest 
that the environmental density dependence of the H2 content is ex- 
tremely subtle except at the lowest densities. We do detect molecu- 
lar gas in galaxies covering the full range of densities and dynamical 
masses probed in our sample. Thus, if the rejuvenation episodes are 
indeed connected to the present-day H2 content, we should expect 
that they are not exclusively limited to low mass galaxies in low 
density environments. A model which fits all of this data should 
also allow for significant rejuvenation rates in high mass galaxies 
and in high density environments. 



6.2 Molecular masses of early-type galaxies 

The cumulative H2 and M(H2)/La- distribution functions in Fig- 
ures [7] and QT] remind us that the majority of the galaxies studied 
here are not detected in CO emission, so that we are really only 
studying the gas-rich ends of the distribution functions. 

Given that the galaxies studied here are early-type galaxies, 
the molecular masses detected are surprisingly high in some cases. 
Three of the sample galaxies have molecular masses greater than 
10 9 Mq (assuming their true CO-H2 conversions are not too differ- 
ent from that used here). These galaxies have comparable amounts 
of m olecular gas to the Milky Way and several times more than in 
M31 dPame et alJI 19931) , yet they are clearly not spirals. The high- 
est molecular masses found in this sa mple are also compa rable to 
that of the merger remnant NGC 7252 dDupraz et al.lll990l) and of 
some lenticular polar ring galaxies dGallettaetal.lll997l) . Nine of 
the sample have molecular masses greater than 5 x 10 8 Mq, which 



is comparabl e to the molecular m ass of the radio galaxy Fornax A 
(NGC 1316; iHorellou et ai]|200ll) . The mass of molecular gas in 
Cen A (NGC 5128) is also a few 10 s Mq dMorgantill2"oiol) . at an 
assumed distance of 4.2 Mpc dTonry et alj200ll) . These two galax- 
ies meet all of the ATLAS 3D selection criteria except for the Dec- 
lination limit, so they could have been in the sample if they were 
elsewhere on the sky. The comparison to Cen A is especially no- 
table because, in spite of its most flamboyant optical appearance, 
Cen A would not have been the most CO-rich galaxy in the sample. 
At M(H2)/Lk ~4xl0~ 3 ,it would also have been solidly in the 
middle of our detected M(H 2 )/L K values (3 x 1CT 4 to 6.3 x 1CT 2 
Mq/Lq). Evidently the molecular gas content of a galaxy does not 
serve as a reliable predictor of its optical morphology, nor does the 
optical morphology reliably predict the molecular content. 

The molecular gas masses detected here are, however, not 
quite in the same league as that of Per seus A (NGC 1275) , 
which has some 10 10 M of cold H 2 dSalome et all l2006h . 
NGC 1275 is the central galaxy of the Perseus Cluster, and 
iLim. Ao. & Dinh-V-Trund d2008l) have recently shown that its 
molecular gas is distributed in filaments which are free-falling to- 
wards the center of the potential. The behaviour is kinematically 
consistent with the molecular gas having cooled and condensed out 
of the galaxy's X-ray halo. But this behaviour is probably not taking 
place in the bulk of the CO-detected ATLAS 3D galaxies, since the 
preponderance of double-peaked CO profiles suggests relaxed discs 
in regular rotation. If the molecular gas in the ATLAS 3D galaxies 
did originate through condensation, it must have happened at least 
a crossing time ago, so that the gas has settled. 

Both the CO detection rate and the H2 mass distribution func- 
tion are independent of the stellar luminosity, for early-type galax- 
ies. When considering all types of galaxies, of course, there is a very 
strong dependence of galaxy properties on the stellar mass; galaxies 
with stellar masses above 3 x 10 10 Mq tend to be old spheroids, 
hence poor in cold gas, and l ess massive galaxies te nd to be HI- 
and CO-rich young discs (e.g. Kauffmann et al. 2003). Among the 
red sequence galaxies, Ivan der Wei et al] d2009h argue that a stel- 
lar mass of 10 Mq is a watershed above which there are very 
few flattened, disc-like systems. However, among the red sequence 
galaxies, we find no evidence that a stellar mass of 3x 10 10 Mq oi- 
IO 11 Mq has any particular significance for a galaxy's molecular 
content. 



6.3 On the origin of the molecular gas in early-type galaxies 

Our results show that at least 22% of all early-type galaxies contain 
molecular gas, and Figure [TT] presents the inferred M(H2)/Lr- dis- 
tribution for M(H 2 )/La- > 10~ 4 Mq/Lq . It is worthwhile to recall 
that this volume-limited survey, complete down to Mk = —21.5, 
has no selection on either FIR luminosity or B magnitude, so that 
it finally eliminates selection bias towards the galaxies which are 

more likely to have star formation activity. 

Using interferometric HI observations, lOosterloo et"aL1 d2010h 
have recently argued that ongoing accretion of cold atomic gas is 
common for field early-type galaxies. In some cases the atomic gas 
retains a large angular momentum and forms structures which are 
tens of kpc in extent, but in other cases the atomic gas is found 
in th e central few kpc of th e galaxy, coincident with a molecular 
disc. lOosterloo et al. I d2010l) thus propose that at least some of the 
molecular gas in early-type galaxies may have been accreted in this 
manner, possibly in the form of atomic gas, settling into a regu- 
lar disc and converting into molecular gas as its density increases. 
Future comparisons of molecular, atomic, ionized gas and stellar 
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kinematics in the ATLAS galaxies should help quantify the role 
of cold accretion as an origin for the molecular gas. 

There is a notable decrease in the CO detection rate of the slow 
or "non-regular" rotators. Since this effect remains even after con- 
trolling for the mass distributions of slow and fast rotators, it must 
be driven by Xr rather than by mass. In other words, the molecular 
content is more strongly correlated with the types of orbits the stars 
occupy than with the number of stars. 

The rate of internal stellar mass loss in a galaxy should natu- 
rally be independent of the angular momentum of the stars. How- 
ever, the mass loss material should retain some memory of the spe- 
cific angular momentum of its progenitor stars. If the material suf- 
fers shock heating, cooling, and condensation to form molecular 
gas, it may lose some specific angular momentum in the process, but 
we would generally expect the mass loss material in a slowly rotat- 
ing galaxy to have smaller specific angular momentum than that in 
a fast rotator. If this scenario accounts for the bulk of the molecular 
gas in early-type galaxies, then, we might expect the molecular gas 
discs in the galaxies of small (intrinsic) A_r to be more compact than 
in the large-Afl galaxies. Indeed, the molecular gas in the galaxies 
of small A_r might have dropped to the nucleus and been consumed, 
instead of forming a kpc-scale disc. From a different perspective, if 
the formation processes of slow rotators are different from those of 
fast rotators, it might also be the case that the formation of a slow 
rotator inevitably destroys any cold gas disc. In short, there are at 
least two scenarios which can explain a lack of cold gas in the slow 
rotators. 

The detection of CO in Virgo Cluster early-type galaxies pro- 
vides a useful perspective on the issue of the origin of the molec- 
ular gas. We have argued above that the CO detection rate among 
ATLAS 3D Virgo Cluster members may be modestly lower than the 
detection rate outside the cluster, but the statistical significance of 
the result is not great. Furthermore, the M(H2) and M(H2)/Lx dis- 
tributions inside and outside the cluster are consistent. We have also 
argued that the CO-detected galaxies are virialized in the cluster po- 
tential and deeply embedded within the hot intracluster gas. Despite 
a long residence in the cluster, their molecular content is not much 
lower than in early-type galaxies outside the cluster. 

Both the internal and external scenarios for the origin of 
molecular gas seem significantly more problematic when galaxies 
are already deep in the cluster potential. Cold mode accretion of gas 
through filaments may be viable for relatively isolated galaxies, but 
gas falling into the cluster should more likely end up as hot gas in 
the general cluster potential than as cold gas in a specific galaxy. 
Acquisition of molecular gas through major or minor mergers also 
should be rare, both because the cluster's high velocity dispersion 
decreases t he merger rate and because mergers tend to fling gas to 
large radii (Barnes 2002) where it would be highly vulnerable to 
stripping. These general expectation s are supported by the N-body 
plus hydrodynamical simulations of iTonnesen et alj d2007b . which 
find that galaxies commonly accrete gas when they are far in the 
outskirts of a cluster but rarely accrete gas after they pass through 
the cluster's virial radius. In addition, the internal stellar mass loss 
scenario also seems difficult because that material is expected to 
shock heat to X-ray temperatures (e.g. iMathews & Brighentill2003l) 
and the concomitant low densities will again make the material vul- 
nerable to stripping. More detailed simulations of the evolution of 
hot gas in cluster members (not cluster-dominant or cluster-centered 
galaxies) would help to quantify this impression. In short, these 
considerations suggest that the CO-detected Virgo Cluster early- 
type galaxies probably did not acquire their cold gas after enter- 
ing the cluster. This in turn suggests that they have retained their 



molecular gas since they fell into the cluster, and that the gas has 
not cycled through a low density phase in that time. 

6.4 On the removal of H2 

If the CO-rich Virgo Cluster early-type galaxies are indeed virial- 
ized, so that they have been in the cluster for at least a relaxation 
time yet have still retained their molecular gas, we conclude that it 
must be difficult to remove the molecular gas entirely. For a cluster 
such as Virgo, with several hundreds to thous ands of galaxies, the 
relaxation time is at least 10 crossing times dBinnev & Tremaing 
2008) and the crossing times are on the order of a Gyr. Thus we 
infer that the virialized, CO-rich early-type galaxies have probably 
retained their gas over at least several Gyr in the cluster. 

It has been suspected for many years that it would be more 
difficult to remove the molecular gas from cluster galaxies than to 
remove the atomic gas. The reasoning is straightforward: the molec- 
ular gas has much higher volume and surface densities than atomic 
gas, and in addition it tends to be deeper in the galaxy's potential 
well, both of which make it more tightly bound. Individual galaxies 
do somet imes show evidence that their molecular hydrogen is bein g 
stripped Jsivanandam. Rieke & Riekd l20ld : IVollmer et al.ll2008h . 
but in general the observations of molecular gas in Virgo Cluster 
spirals have suggested that H2 deficiencies are su btle, even when 
HI deficiencies are strong dKennev & Y oung 1989). The same now 
appears to be true for Virgo Cluster early-type galaxies, where H2 
deficie ncies are subtle to nonexistent whereas HI deficien cies are 
strong (di Serego Alighie ri et al. I l2007t iGrossi et al.l l2009). In this 
context, the Virgo Cluster early-type galaxies provide a useful per- 
spective because (as we argued above) they have been residing in 
the cluster longer than the spirals. They therefore provide longer 
time-baseline views into the gas removal processes than spirals do. 

Recent adap t ive-me sh hydrodynamical simulations of 
ITonnesen & Brvanl d2009l) suggest that when a galaxy first falls 
into the intracluster medium, it suffers a brief but intense period 
of stripping when most of the loosely-bound material is removed. 
In the simulations, dense molecular gas can be retained through 
this period, in agreement with the implicati ons of our work. 
Quantitatively, the standard fGunn &Gottl ( ll972h stripping criterion 
requires a ram pressure ~ V^'Sgas/R- For typical values we take 
a circular velocity V c of 200 km s -1 , molecular gas surface density 
Egas at 100 M© pc -2 , and a ra dius R of 1 kpc, bas ed on the 
interferometric CO observations of ICrocker et al.l ( 2010l). The ram 
pressure needed is then ~ 3 x 10~ 9 dyne cm' - 2 . IVollmej|2009l) 
finds that these kinds of pressures require a galaxy to approach 
within ~ 50 kpc of the Virgo Cluster center even if the velocity at 
closest approach is 2000 km s~ . Thus the molecular gas should be 
retained except on very h igh eccentricity orbits within the cluster. 
ITonnesen & Brvanl (|2009j) also mention that angular momentum 
transfer to the non-rotating cluster gas can gradually drive the 
remaining molecular gas towards the center of its host galaxy. For 
all of these reasons, the CO-rich Virgo Cluster early-type galaxies 
could motivate longer simulation runs in an attempt to test whether 
these processes can transform an exponential gas disc (typical 
of spirals) into a rather compact, sharp-edged molecular disc 
such as those observed in NGC4459, NGC4526 and NGC 4477 
dYoung. Bureau & Ca ppellari 20081; ICrocker et al.ll2010l) . 

6.5 On the formation of fast rotators 

While our early-type galaxies are today rather red in colour, they 
must have been blue during the epoch when their stars were form- 
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ing. The mass distributions and morphologies of their progenitors 
are, however, the subjects of much debate. One major question is 
whether or not the fast rotator galaxies passed throug h a stage as 
spiral galaxies. Some authors suggest not. For example. iDekel et alJ 
( 2009) propose that if the matter accreting through cold streams into 
a protogalactic halo is highly clumped, the dynamical interactions 
between the clumps will counteract the tendency of cold gas to set- 
tle into a disc. In this manner a spheroid-dominated galaxy could 
form without ever passing through a disc-dominated stage, even if 
the bulk of the accreted baryons are in the form of cold gas. 

Other popular scenarios propose to form fast rotator early- 
type galaxies by converting spirals. One aspect of the conversion 
must involve quenching the star formation activity (so as to red- 
den and age the stellar population); another aspect is to effect the 
morphological change (dispersal of the spiral structure, disc thick- 
ening and/or bulge growth). It is often assumed that the removal 
of the cold gas will accomplish these objectives. Turning off new 
star formation would remove the continuous input of a dynamically 
cold stellar population, so that dynamical heating and thickening 
would naturally erase the non-axisymmetric spiral structures. Re- 
moval of the cold gas may not be nec essary for this conversion pro- 
cess, though. iBournaud et al.l d2007l) have suggested that multiple 
minor mergers (with mass ratios of 10: 1 and greater) will gradually 
thicken a spira l's stellar dis c and increase its central concentration 
(Sersic index). Martig et all ( l2009h also remind us that if the gravi- 
tational potential can be sufficiently deepened, a cold gas disc could 
be simply rendered incapable of star formation activity. Our current 
work suggests that, indeed, some cluster galaxies can retain sig- 
nificant amounts of molecular gas even after being resident in the 
cluster for several Gyr. 

We cannot, at present, confirm or refute the speculation that 
fast rotator early-type galaxies in the Virgo Cluster are the stripped, 
heated, and/or quenched remnants of spirals. We simply comment 
that if this is the case, the paradigm should not necessarily require 
all of the molecular gas to be removed in the conversion process. 
In a purely hydrodynamic stripping scenario, it would be interest- 
ing to test whether the star formation quenching and morphologi- 
cal change can be accomplished quickly enough by removing only 
some of the molecular gas and simply rearranging the rest. In the 
merger, morphological quenching, and clumpy accretion scenarios, 
it would be valuable to have statistical predictions for the proper- 
ties (masses and spatial extents) of any cold gas discs which might 
remain. 



7 SUMMARY 

We present a 12 CO J=l-0 and J=2-l search in the galaxies of the 
ATLAS 3D sample, a complete volume-limited sample of early-type 
galaxies with —21.5 ^ Mk ^ — 26 and with morphologies veri- 
fied by inspection of optical (mostly SDSS) images. CO data from 
the IRAM 30m telescope are collected for 259 of the 260 mem- 
bers, including 204 new observations and the remainder collected 
from the literature. The 3a upper limit for a sum over a 300 km s _1 
linewidth corresponds to an H2 mass ~ 1 x 10 7 Mq for the nearby 
sample members (11 Mpc) and a mass ~lx 10 8 Mq for the most 
distant members (40 Mpc). A few targets are observed significantly 
more deeply than that. 

The detected CO line intensities correspond to H2 masses of 
logM(H 2 )/M Q = 7.10 to 9.29 and M(H 2 )/L K = 2.8 x 10" 4 to 
0.063 M /L . The CO detection rate is 56/259 = 0.22±0.03, a rate 
known to be a lower limit since (1) molecular gas outside the cen- 



tral 22" has not been observed, (2) we expect to be incomplete for 
linewidths greater than 300 km s" 1 , and (3) compact distributions 
suffer strong beam dilution. Many of the detections are brighter in 
CO(2-l) than in CO(1-0). The line profiles are often double-horned 
and/or asymmetric. These observations collectively suggest that the 
molecular gas is often found in regular rotating discs. It is some- 
times in compact structures (smaller than the 12" CO(2-l) beam, 
which is 1.1 kpc at 20 Mpc) and sometimes in structures com- 
parable to or larger than the CO(1-0) beam (2.2 kpc at 20 Mpc). 
The linewidth distribution peaks at smaller widths than expected, 
suggesting that we have some cases in which the molecular gas 
reaches the galaxy's asymptotic circular velocity and some in which 
it doesn't. Molecular surface densities, averaged over the inner 2 to 
4 kpc of the galaxies, are comparable to those found in spirals. Inter- 
ferometric observations will be necessary in order to quantify those 
general impressions about the typical gas distributions. 

There is a strong correlation between the presence of molec- 
ular gas and the presence of dust, blue features, and young stellar 
ages seen in H/3 absorption. Thus, the molecular gas that we detect 
is often engaged in star formation. 

The CO detection rate is not a function of luminosity over the 
range —21.5 > Mk ^ —26, and the H2 mass distribution function 
is also independent of stellar luminosity. The CO detection rate is 
a strong function of the global velocity width a e , with galaxies of 
<7 C < 100 km s _1 more than twice as likely to be detected. How- 
ever, the dependence on the dynamical mass inferred from mod- 
elling the stellar kinematics is considerably more subtle. These con- 
siderations suggest that the a e effect is probably a combination of 
a downsizing trend and an observational bias in favor of detecting 
CO in face-on galaxies. 

The CO detection rate is strongly correlated with the stellar 
specific angular momentum and internal kinematic substructure, 
such that slow rotators and non-regular rotators have much less 
molecular gas than their fast and regular rotating counterparts. This 
result remains true after controlling for the a e and dynamical mass 
dependences noted above. We suggest that either the assembly pro- 
cesses which built the slow rotators also destroyed or removed their 
molecular gas, or they are less likely to re-acquire molecular gas 
after assembly. 

Approximately half of the fast rotators at distances less than 
24 Mpc are members of the Virgo Cluster and half are not. The CO 
detection rate is modestly lower in cluster members than in non- 
members; however, it is important to note that the drop in detection 
rate is only a ~ la effect. There is no measurable difference in the 
H2 or M(H2)/Lx distribution functions of cluster and non-cluster 
galaxies within 24 Mpc. CO-rich Virgo Cluster early- type galaxies 
are virialized in the cluster potential and they do not avoid the center 
of the cluster (unlike cluster spirals). Apparently they have retained 
their molecular gas through extended residences in the cluster. 

The rarest, most gas-rich galaxies (M(H 2 )/Lk> 0.02 
Mq/Lq) are found in relatively low-density environments at dis- 
tances beyond 24 Mpc. These results suggest that the CO-detected 
Virgo Cluster galaxies have CO-detected analogs of the same 
M(H 2 )/La- in the field, but there is in addition a small population of 
CO-rich objects in the field that is not found in the Virgo Cluster. 
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Table 4. CO intensities and H2 masses. 



Galaxy 


rms (1-0) 


rms (2-1) 


Vel. range 


1(1-0) 


1(2-1) 


Ratio 


Source log M(H 2 ) 




mK 


mK 


km s _1 


Kkms" 1 


Kkms" 1 


(2-l)/(l-0) 


M 



IC0560 


3.86 


9.40 


1703.0 


2003.0 


0.106(0.347) 


-0.830 ( 0.758) 






< 7.67 


IC0598 


4.19 


5.72 


2106.0 


2406.0 


-0.007 ( 0.457) 


-0.340 ( 0.525) 




j 


< 8.02 


IC0676 


4.84 


4.92 


1310.0 


1516.0 


11.641 ( 0.429) 


16.654 ( 0.343) 


1.43 (0.06) 




8.63 (0.02) 


IC0719 


2.80 


2.74 


1619.0 


2022.0 


3.472(0.361) 


5.798 ( 0.311) 


1.67 (0.20) 


j 


8.26 (0.04) 


IC0782 


3.05 


3.86 


2274.0 


2574.0 


0.337 ( 0.343) 


0.072 ( 0.333) 




j 


< 7.92 


IC1024 


4.75 


6.35 


1359.0 


1619.0 


11.419 ( 0.428) 


15.631 (0.401) 


1.37 (0.06) 


j 


8.61 (0.02) 


IC3631 


3.23 


3.17 


2672.0 


2972.0 


0.213 (0.268) 


-0.006 ( 0.235) 




1 


< 7.94 


NGC0448 


3.72 


4.57 


1758.0 


2058.0 


-0.399 ( 0.345) 


-0.767 ( 0.395) 




1 


< 7.74 


NGC0474 


















< 7.68 


NGC0502 


3.08 


3.57 


2374.0 


2674.0 


-0.344 ( 0.320) 


0.110(0.291) 




1 


< 7.88 
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3.50 


4.27 


2232.0 


2323.0 


0.477 ( 0.155) 


0.348 (0.165) 


< 1.04 


! 


7.48 (0.12) 


NGC0516 
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0.544 ( 0.442) 
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< 7.82 


NGC0524 


3.07 
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2618.0 


2.843 ( 0.348) 
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2289.0 


0.017 ( 0.330) 


0.658 ( 0.367) 
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< 7.75 
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< 7.52 
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< 7.47 
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4.41 


6.13 


2253.0 


2603.0 


17.354 ( 0.524) 


28.760 ( 0.560) 


1.66 (0.06) 


j 


9.07 (0.01) 


NGC1248 


3.11 


3.55 


2067.0 


2367.0 


0.307 ( 0.287) 


0.055 (0.281) 




j 


< 7.68 


NGC1266 


5.82 


4.33 


1750.0 


2500.0 


34.959 ( 1.068) 


104.397 ( 0.926) 


2.99 (0.09) 


j 


9.28 (0.01) 


NGC1289 


3.41 


4.49 


2642.0 


2942.0 


0.263 ( 0.285) 


-0.003 ( 0.343) 




j 


< 7.89 


NGC1665 


3.62 


5.53 


2595.0 


2895.0 


0.506 ( 0.349) 


0.112 ( 0.394) 




j 


< 7.95 


NGC2481 


3.36 


5.79 


2007.0 


2307.0 


0.030 ( 0.326) 


0.136 ( 0.492) 




j 


< 7.79 


NGC2549 


















< 7.06 


NGC2577 


3.51 


5.45 


1912.0 


2212.0 


-0.113 (0.294) 


-0.350 ( 0.377) 




j 


< 7.71 


NGC2592 


3.55 


4.10 


1829.0 


2129.0 


0.319 ( 0.305) 


0.905 ( 0.348) 




j 


< 7.54 


NGC2594 


3.83 


3.18 


2212.0 


2512.0 


0.294 ( 0.302) 


0.046 ( 0.258) 




! 


< 7.83 


NGC2679 


















< 7.87 


NGC2685 


2.78 


7.20 


746.0 


931.0 


1.158 ( 0.225) 


0.679 ( 0.425) 


< 1.10 


! 


7.29 (0.08) 


NGC2695 


















< 8.01 


NGC2698 


2.94 


0.00 


1750.0 


2050.0 


-0.003 ( 0.237) 


0.000 ( 0.000) 




j 


< 7.50 


NGC2699 


















< 7.54 


NGC2764 


5.21 


7.18 


2514.0 


2944.0 


16.198 (0.580) 


24.398 ( 0.847) 


1.51 (0.08) 


j 


9.19 (0.02) 


NGC2768 


















7.64 (0.07) 


NGC2778 


3.66 


5.21 


1875.0 


2175.0 


0.183 (0.330) 


-0.719(0.412) 




j 


< 7.48 


NGC2824 


1.95 


5.23 


2552.0 


2996.0 


4.460 ( 0.334) 


6.175 (0.528) 


1.38 (0.16) 




8.65 (0.03) 


NGC2852 


4.15 


4.37 


1631.0 


1931.0 


0.176(0.325) 


0.615 (0.359) 






< 7.68 


NGC2859 


3.71 


4.73 


1540.0 


1840.0 


0.319 (0.309) 


0.220 ( 0.364) 






< 7.61 


NGC2880 


3.43 


7.69 


1404.0 


1704.0 


-0.088 ( 0.336) 


-0.271 (0.511) 






< 7.44 


NGC2950 


3.54 


4.89 


1172.0 


1472.0 


-0.115 (0.347) 


-0.020 ( 0.465) 






< 7.12 


NGC2962 


3.73 


8.97 


1817.0 


2117.0 


0.074 ( 0.336) 


-0.309 (0.712) 






< 7.85 


NGC2974 
















2 


< 7.65 


NGC3032 
















2 


8.41 (0.01) 


NGC3073 
















4 


7.52 (0.07) 


NGC3098 


3.59 


4.03 


1247.0 


1547.0 


0.854 ( 0.305) 


0.222 ( 0.341) 




1 


< 7.47 


NGC3156 
















2 


7.67 (0.09) 
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Table 4 (cont'd) 



Galaxy 


rms (1-0) 


rms (2-1) 


Vel. range 


^(1-0) 


^(2-1) 


Ratio 


Source 


log M(H 2 ) 




mK 


mK 


km s 




Kkms" 1 


Kkms" 1 


(2-l)/(l-0) 




M 


NGC3182 


3.52 


2.81 


1942.0 


2239.0 


3.007 ( 0.346) 


4.166 (0.265) 


1.39 (0.18) 


1 


8.33 (0.05) 


NGC3193 


4.04 


4.19 


1231.0 


1531.0 


-0.081 (0.405) 


-0.383 (0.315) 




1 


< 7.91 


NGC3226 


3.52 


3.72 


1165.0 


1465.0 


0.306(0.271) 


0.635 ( 0.294) 




1 


< 7.41 


NGC3230 


3.52 


7.08 


2645.0 


2945.0 


-0.374 ( 0.328) 


-1.105 (0.571) 




1 


< 8.00 


NGC3245 


3.67 


4.76 


1320.0 


1495.0 


0.736 ( 0.228) 


0.125 (0.256) 


< 1.04 


1 


7.27 (0.12) 


NGC3248 


3.47 


5.34 


1331.0 


1631.0 


0.095 ( 0.321) 


-0.672 ( 0.483) 




1 


< 7.55 


NGC3301 


3.57 


3.13 


1189.0 


1489.0 


0.526(0.301) 


0.163 (0.271) 




1 


< 7.46 


NGC3377 
















2 


< 6.96 


NGC3379 
















2 


< 6.72 


NGC3384 
















2 


< 7.11 


NGC3400 


3.94 


3.90 


1291.0 


1591.0 


-0.363 ( 0.384) 


-0.090 ( 0.309) 




1 


< 7.63 


NGC3412 


4.41 


3.36 


710.0 


1010.0 


-0.238 (0.413) 


-0.283 ( 0.269) 




1 


< 6.96 


NGC3414 
















2 


< 7.19 


NGC3457 


3.41 


4.58 


998.0 


1298.0 


0.489 ( 0.303) 


-0.070 ( 0.306) 




1 


< 7.35 


NGC3458 


4.06 


3.04 


1727.0 


2027.0 


-0.674 ( 0.307) 


0.215 (0.223) 




1 


< 7.73 


NGC3489 
















2 


7.20 (0.06) 


NGC3499 


3.57 


3.21 


1385.0 


1685.0 


0.861 ( 0.325) 


0.220 ( 0.249) 




1 


< 7.62 


NGC3522 
















7 


< 7.28 


NGC3530 


4.14 


3.99 


1744.0 


2044.0 


0.724 ( 0.341) 


-0.398 ( 0.305) 




1 


< 7.78 


NGC3595 


3.24 


2.67 


2027.0 


2327.0 


0.873 (0.311) 


-0.043 ( 0.235) 




1 


< 7.84 


NGC3599 


2.70 


3.07 


687.0 


893.0 


0.638 ( 0.221) 


0.959 (0.183) 


> 1.44 


1 


7.36 (0.08) 


NGC3605 
















4 


< 7.48 


NGC3607 
















3 


8.42 (0.05) 


NGC3608 
















2 


< 7.58 


NGC3610 


3.46 


5.04 


1557.0 


1857.0 


-0.103 (0.320) 


0.026 ( 0.470) 




1 


< 7.40 


NGC3613 


3.37 


3.98 


1901.0 


2201.0 


-0.356 ( 0.309) 


0.312 (0.293) 




1 


< 7.66 


NGC3619 


3.84 


4.12 


1303.0 


1885.0 


4.318(0.510) 


3.806 ( 0.677) 


0.88 (0.19) 


1 


8.28 (0.05) 


NGC3626 


5.25 


4.25 


1202.0 


1731.0 


6.998 ( 0.658) 


11.285 (0.569) 


1.61 (0.17) 


1 


8.21 (0.04) 


NGC3630 


3.67 


4.06 


1349.0 


1649.0 


-0.476 ( 0.350) 


0.482 ( 0.345) 




1 


< 7.60 


NGC3640 


2.85 


3.24 


1148.0 


1448.0 


0.291 ( 0.306) 


0.382 ( 0.265) 




1 


< 7.59 


NGC3641 


3.72 


4.77 


1630.0 


1930.0 


0.729 ( 0.373) 


0.036 ( 0.378) 




1 


< 7.66 


NGC3648 


3.51 


5.38 


1820.0 


2120.0 


-0.135 (0.319) 


0.008 ( 0.372) 




1 


< 7.77 


NGC3658 


3.93 


4.14 


1889.0 


2189.0 


0.233 ( 0.339) 


0.582 ( 0.360) 




1 


< 7.82 


NGC3665 


3.48 


4.41 


1737.0 


2432.0 


12.159 (0.656) 


14.247 ( 0.836) 


1.17(0.09) 


1 


8.91 (0.02) 


NGC3674 


3.70 


3.73 


1905.0 


2205.0 


0.125 (0.298) 


0.272 ( 0.258) 




1 


< 7.78 


NGC3694 


3.63 


5.17 


2093.0 


2393.0 


0.340 ( 0.360) 


0.488 ( 0.435) 




1 


< 7.91 


NGC3757 


3.88 


3.02 


1095.0 


1395.0 


-0.177 (0.326) 


-0.686 ( 0.248) 




1 


< 7.48 


NGC3796 


2.49 


3.16 


1110.0 


1410.0 


-0.150 (0.340) 


-0.001 ( 0.270) 




1 


< 7.51 


NGC3838 


4.32 


3.45 


1158.0 


1458.0 


-0.547 ( 0.338) 


-0.016 ( 0.272) 




1 


< 7.53 


NGC3941 


3.17 


3.55 


780.0 


1080.0 


0.090 ( 0.300) 


-0.010 ( 0.272) 




1 


< 6.89 


NGC3945 


4.07 


2.98 


1131.0 


1431.0 


-0.122(0.320) 


0.046 ( 0.241) 




1 


< 7.50 


NGC3998 


3.62 


3.57 


898.0 


1198.0 


0.053 ( 0.338) 


-0.033 ( 0.270) 




1 


< 7.06 


NGC4026 


3.03 


2.82 


835.0 


1135.0 


-0.169(0.310) 


0.015 (0.208) 




1 


< 6.99 


NGC4036 


2.79 


2.29 


1169.0 


1670.0 


0.372 ( 0.492) 


3.676 ( 0.323) 


> 2.49 




8.13 (0.04) 


NGC4078 


4.51 


7.41 


2396.0 


2696.0 


-0.088 ( 0.360) 


-1.343 (0.573) 






< 7.98 


NGC4111 


















7.22 (0.09) 


NGC4119 


4.22 


7.34 


1543.0 


1767.0 


4.537 ( 0.320) 


5.787 (0.435) 


1.28 (0.13) 




7.88 (0.03) 


NGC4143 


















< 7.20 


NGC4150 


















7.82 (0.03) 


NGC4168 


3.61 


6.51 


2136.0 


2436.0 


0.217 (0.314) 


-0.219(0.514) 






< 7.74 
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Table 4 (cont'd) 



Galaxy 


rms (1-0) 


rms (2-1) 


Vel. range 


^(1-0) 


^(2-1) 


Ratio 


Source 


log M(H 2 ) 




mK 


mK 


km s 




Kkm s" 1 


Kkms" 1 


(2-l)/(l-0) 




M 


NGC4179 


3.95 


5.57 


1150.0 


1450.0 


0.519 (0.383) 


0.534 ( 0.388) 




1 


<7.28 


NGC4191 


3.40 


4.11 


2496.0 


2796.0 


0.493 (0.310) 


-0.005 ( 0.263) 




1 


<7.94 


NGC4203 
















3 


7.39 (0.05) 


NGC4215 


3.85 


4.36 


1861.0 


2161.0 


-0.007 ( 0.370) 


-0.451 (0.385) 




1 


<7.83 


NGC4233 


4.16 


9.28 


2156.0 


2456.0 


0.109 (0.367) 


0.417 ( 0.734) 




1 


<7.89 


NGC4249 


3.84 


3.47 


2468.0 


2768.0 


0.264 ( 0.340) 


0.168 (0.281) 




1 


<7.97 


NGC4251 
















3 


< 7.1 1 


NGC4255 


3.63 


3.57 


1845.0 


2145.0 


0.220 ( 0.340) 


0.250 ( 0.252) 




1 


<7.78 


NGC4259 


4.11 


7.04 


2347.0 


2647.0 


0.523 ( 0.368) 


0.996 ( 0.589) 




1 


<7.97 


NGC4261 
















2 


<7.68 


NGC4262 
















2 


<7.07 


NGC4264 


3.54 


3.42 


2368.0 


2668.0 


0.194 (0.342) 


-0.024 ( 0.272) 




1 


<7.94 


NGC4267 


4.18 


4.80 


871.0 


1171.0 


0.198 (0.314) 


-0.450 ( 0.432) 




1 


< 7.16 


NGC4268 


3.72 


4.31 


1884.0 


2184.0 


0.140(0.372) 


-0.223 ( 0.347) 




1 


<7.83 


NGC4270 
















2 


<7.79 


NGC4278 
















7 


<7.45 


NGC4281 


3.12 


3.78 


2302.0 


3073.0 


1.286 (0.693) 


1.572(0.602) 




1 


<7.88 


NGC4283 
















4 


7.10(0.09) 


NGC4324 


3.57 


3.71 


1536.0 


1852.0 


2.955 ( 0.334) 


1.977 (0.290) 


0.67 (0.12) 


1 


7.69 (0.05) 


NGC4339 


3.70 


3.82 


1116.0 


1416.0 


0.238 ( 0.305) 


-0.632 ( 0.286) 




1 


< 7.15 


NGC4340 


4.17 


4.74 


783.0 


1083.0 


-0.004 ( 0.349) 


-0.466(0.371) 




1 


<7.33 


NGC4342 


3.93 


3.67 


611.0 


911.0 


0.121 (0.346) 


0.326 ( 0.297) 




1 


<7.24 


NGC4346 


3.79 


4.11 


682.0 


982.0 


0.037 ( 0.374) 


-0.534 ( 0.308) 




1 


<7.12 


NGC4350 


3.77 


3.85 


1060.0 


1360.0 


0.350 ( 0.352) 


0.501 (0.349) 




1 


< 7.18 


NGC4365 
















2 


<7.62 


NGC4371 


3.58 


6.85 


783.0 


1083.0 


-0.396 ( 0.370) 


-0.711 (0.634) 




1 


<7.29 


NGC4374 
















2 


<7.23 


NGC4377 


3.95 


3.78 


1193.0 


1493.0 


0.424(0.315) 


-0.458 ( 0.276) 




1 


<7.26 


NGC4379 


3.64 


6.74 


924.0 


1224.0 


-0.455 ( 0.338) 


0.071 (0.503) 




1 


< 7.19 


NGC4382 
















2 


<7.39 


NGC4387 
















2 


<7.39 


NGC4406 
















5 


<7.40 


NGC4417 


3.14 


3.45 


678.0 


978.0 


-0.357 ( 0.359) 


-0.681 (0.306) 




1 


< 7.22 


NGC4425 


2.87 


3.16 


1758.0 


2058.0 


0.027 (0.318) 


0.562 ( 0.294) 




1 


<7.20 


NGC4429 


3.03 


5.14 


808.0 


1405.0 


6.769 ( 0.497) 


7.334 ( 0.747) 


1.08 (0.14) 


1 


8.05 (0.03) 


NGC4434 


3.77 


5.20 


920.0 


1220.0 


-0.119(0.431) 


-0.045 ( 0.526) 




1 


<7.60 


NGC4435 


2.47 


5.36 


569.0 


1022.0 


4.361 (0.389) 


5.667 ( 0.686) 


1.30 (0.20) 


1 


7.87 (0.04) 


NGC4442 


3.84 


7.72 


397.0 


697.0 


0.297 (0.311) 


0.678 ( 0.574) 




1 


<7.12 


NGC4452 


3.94 


3.45 


38.0 


338.0 


0.020 ( 0.346) 


0.364 ( 0.262) 




1 


< 7.19 


NGC4458 
















2 


< 7.31 


NGC4459 
















2 


8.24 (0.02) 


NGC4461 


3.72 


4.74 


1774.0 


2074.0 


-0.256 (0.319) 


0.333 ( 0.399) 




1 


<7.20 


NGC4472 


4.12 


5.03 


831.0 


1131.0 


-0.472 ( 0.330) 


0.050(0.417) 




1 


<7.25 


NGC4473 
















2 


<7.07 


NGC4474 


3.47 


7.41 


1461.0 


1761.0 


0.359 ( 0.329) 


0.365 ( 0.482) 




1 


< 7.16 


NGC4476 
















6 


8.05 (0.04) 


NGC4477 
















2 


7.54 (0.06) 


NGC4478 


2.91 


4.31 


1225.0 


1525.0 


-0.379 ( 0.363) 


-0.109 (0.377) 




1 


<7.28 


NGC4483 


3.82 


3.17 


756.0 


1056.0 


0.644(0.314) 


-0.153 (0.226) 




1 


< 7.20 


NGC4486 
















2 


< 7.17 
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Table 4 (cont'd) 



Galaxy 


rms (1-0) 


rms (2-1) 


Vel. range 


1(1-0) 


1(2-1) 


Ratio 


Source 


log M(H 2 ) 




mK 


mK 


km s 


-i 


Kkms" 1 


Kkmr 1 


(2-l)/(l-0) 




M 


NGC4486A 



















< 0.00 


NGC4489 


4.03 


5.25 


811.0 


1111.0 


0.215 ( 0.324) 


0.708 (0.431) 




1 


< 7.15 


NGC4494 


3.51 


4.81 


1192.0 


1492.0 


-0.133 (0.349) 


-0.937 ( 0.368) 




1 


< 7.25 


NGC4503 


4.00 


3.40 


1184.0 


1484.0 


-0.083 ( 0.333) 


0.092 ( 0.257) 




1 


< 7.22 


NGC4521 


3.79 


4.44 


2361.0 


2661.0 


-0.030 ( 0.321) 


-0.081 ( 0.342) 




1 


< 7.97 


NGC4526 


















8.59 (0.01) 


NGC4528 


3.27 


4.04 


1228.0 


1528.0 


-0.236 (0.311) 


0.420 ( 0.292) 




1 


< 7.15 


NGC4546 


















< 6.97 


NGC4550 


4.22 


5.16 


309.0 


609.0 


0.234 ( 0.396) 


0.667 (0.571) 




1 


< 7.24 


NGC4551 


3.66 


6.17 


1026.0 


1326.0 


0.100(0.371) 


-0.289 ( 0.495) 




1 


< 7.24 


NGC4552 


















< 7.28 


NGC4564 


4.00 


6.58 


1005.0 


1305.0 


-0.020 ( 0.387) 


0.468 ( 0.497) 




1 


< 7.25 


NGC4570 


















< 7.47 


NGC4578 


3.27 


6.63 


2142.0 


2442.0 


0.384 ( 0.323) 


-0.202 ( 0.558) 




1 


< 7.20 


NGC4596 


3.36 


4.78 


1715.0 


2057.0 


1.246(0.294) 


1.667 (0.390) 


1.34 (0.44) 


1 


7.31 (0.09) 


NGC4608 


3.38 


4.37 


1700.0 


2000.0 


0.039 ( 0.401) 


-0.078 ( 0.359) 




1 


< 7.30 


NGC4612 


3.13 


4.19 


1625.0 


1925.0 


-0.441 (0.313) 


-0.098 (0.331) 




1 


< 7.20 


NGC4621 


3.94 


2.38 


317.0 


617.0 


0.289 ( 0.336) 


-0.141 (0.200) 




1 


< 7.13 


NGC4623 


2.95 


4.16 


1657.0 


1957.0 


-0.394 ( 0.294) 


-0.230 ( 0.289) 




1 


< 7.21 


NGC4624 


3.46 


7.32 


762.0 


1062.0 


0.242 ( 0.400) 


0.092 ( 0.653) 




1 


< 7.30 


NGC4636 


















< 6.87 


NGC4638 


3.73 


5.29 


1002.0 


1302.0 


-0.352 ( 0.357) 


0.384 ( 0.334) 




1 


< 7.30 


NGC4643 


3.76 


3.29 


1220.0 


1673.0 


0.317 (0.481) 


1.111 (0.356) 


>0.77 


1 


7.27 (0.12) 


NGC4649 


8.94 


8.87 


960.0 


1260.0 


0.249 ( 0.508) 


0.214(0.711) 




1 


< 7.44 


NGC4660 


3.73 


6.11 


937.0 


1237.0 


0.131 (0.381) 


-0.050 ( 0.522) 




1 


< 7.19 


NGC4684 


3.86 


5.34 


1420.0 


1844.0 


1.544(0.439) 


6.188 (0.638) 


4.01 (1.21) 


1 


7.21 (0.11) 


NGC4690 


3.45 


4.06 


2615.0 


2915.0 


-0.724 ( 0.350) 


-0.186 (0.327) 




1 


< 8.01 


NGC4694 


6.02 


4.80 


1083.0 


1264.0 


6.182(0.372) 


6.566 ( 0.301) 


1.06 (0.08) 


1 


8.01 (0.03) 


NGC4697 


3.86 


3.84 


1102.0 


1402.0 


0.466 ( 0.305) 


-0.414 ( 0.269) 




1 


< 6.86 


NGC4710 


8.40 


6.89 


896.0 


1368.0 


31.696 (0.923) 


40.130(0.766) 


1.27 (0.04) 


1 


8.72 (0.01) 


NGC4733 


4.66 


3.96 


775.0 


1075.0 


-0.001 ( 0.499) 


0.580 ( 0.308) 




1 


< 7.28 


NGC4753 


4.51 


3.75 


863.0 


1513.0 


11.169 ( 0.679) 


14.203 ( 0.576) 


1.27 (0.09) 


1 


8.55 (0.03) 


NGC4754 


3.39 


3.64 


1201.0 


1501.0 


-0.595 (0.319) 


-0.177 (0.273) 




1 


< 7.18 


NGC4762 


4.14 


7.64 


836.0 


1136.0 


0.087 ( 0.321) 


0.171 (0.517) 




1 


< 7.48 


NGC4803 


3.80 


3.71 


2495.0 


2795.0 


-0.467 ( 0.337) 


0.018 (0.294) 




1 


< 7.98 


NGC5103 


4.13 


3.44 


1123.0 


1423.0 


-0.043 ( 0.383) 


-0.222 ( 0.285) 




1 


< 7.58 


NGC5173 


3.32 


2.51 


2276.0 


2568.0 


2.134(0.288) 


3.109 (0.228) 


1.46 (0.22) 


1 


8.28 (0.06) 


NGC5198 


















< 7.89 


NGC5273 


2.90 


2.87 


961.0 


1214.0 


1.303 (0.226) 


3.309 ( 0.203) 


2.54 (0.47) 


1 


7.31 (0.07) 


NGC5308 


















< 7.88 


NGC5322 


3.26 


4.92 


1630.0 


1930.0 


0.054 ( 0.341) 


-0.083 ( 0.379) 




1 


< 7.76 


NGC5342 


2.57 


3.48 


2039.0 


2339.0 


-0.131 (0.265) 


-0.567 ( 0.262) 




1 


< 7.79 


NGC5353 


3.37 


5.23 


1910.0 


2633.0 


1.358 (0.581) 


2.409 (0.891) 




1 


< 8.12 


NGC5355 


3.95 


7.37 


2194.0 


2494.0 


0.854 ( 0.343) 


1.017 (0.587) 






< 7.94 


NGC5358 


3.47 


4.84 


2262.0 


2562.0 


-0.308 (0.317) 


-0.052 ( 0.365) 






< 7.92 


NGC5379 


4.82 


7.80 


1683.0 


1911.0 


3.914(0.345) 


3.075 ( 0.610) 


0.79 (0.17) 




8.33 (0.04) 


NGC5422 


4.10 


7.52 


1688.0 


1988.0 


-0.326 ( 0.349) 


-0.162 (0.628) 






< 7.78 


NGC5473 


3.27 


5.99 


1872.0 


2172.0 


0.062 ( 0.349) 


0.591 (0.477) 






< 7.85 


NGC5475 


3.80 


4.35 


1521.0 


1821.0 


0.437 ( 0.353) 


0.237 ( 0.384) 






< 7.72 


NGC5481 


3.02 


5.44 


1839.0 


2139.0 


0.133 (0.325) 


0.485 (0.451) 






< 7.60 
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Table 4 (cont'd) 



Galaxy 


rms (1-0) 


rms (2-1) 


Vel. range 


^(1-0) 


^(2-1) 


Ratio 


Source 


log M(H 2 ) 




mK 


mK 


km s 


-i 


Kkms" 1 


Kkms" 1 


(2-l)/(l-0) 




M 


NGC5485 


3.07 


3.99 


1777.0 


2077.0 


0.121 (0.342) 


0.300 (0.318) 




1 


< 7.60 


NGC5493 


3.62 


4.33 


2515.0 


2815.0 


-0.436 ( 0.347) 


0.669 ( 0.369) 




1 


< 7.98 


NGC5500 


3.67 


4.82 


1764.0 


2064.0 


-0.696 ( 0.363) 


-0.019 ( 0.428) 




1 


< 7.82 


NGC5507 


4.23 


2.50 


1701.0 


2001.0 


-0.129 (0.334) 


0.335 ( 0.247) 




1 


< 7.70 


NGC5557 


3.42 


2.54 


3069.0 


3369.0 


-0.566 ( 0.300) 


-0.078 ( 0.238) 




1 


< 7.92 


NGC5574 


3.79 


4.06 


1439.0 


1739.0 


0.016 (0.327) 


-0.392 ( 0.288) 




1 


< 7.51 


NGC5576 


3.06 


3.57 


1356.0 


1656.0 


0.081 (0.354) 


-0.168 (0.272) 




1 


< 7.60 


NGC5582 


3.99 


3.75 


1280.0 


1580.0 


0.499 ( 0.336) 


-0.052 ( 0.274) 




1 


< 7.67 


NGC5611 


3.78 


4.72 


1818.0 


2118.0 


0.675 ( 0.338) 


-0.141 (0.353) 




1 


< 7.57 


NGC5631 


3.13 


4.51 


1794.0 


2094.0 


0.166 (0.357) 


0.744 ( 0.341) 




1 


< 7.68 


NGC5638 


3.08 


3.27 


1502.0 


1802.0 


-0.019(0.331) 


0.026(0.251) 




1 


< 7.60 


NGC5687 


3.71 


2.96 


1993.0 


2293.0 


0.454 ( 0.323) 


0.419(0.261) 




1 


< 7.64 


NGC5770 


3.73 


3.80 


1321.0 


1621.0 


-0.250 ( 0.347) 


-0.250 ( 0.305) 




1 


< 7.34 


NGC5813 
















2 


< 7.69 


NGC5831 
















2 


< 7.85 


NGC5838 
















2 


< 7.56 


NGC5839 


3.94 


2.92 


1070.0 


1370.0 


0.382 ( 0.269) 


0.246 ( 0.217) 




1 


< 7.38 


NGC5845 
















2 


< 7.50 


NGC5846 
















2 


< 7.78 


NGC5854 


3.53 


2.67 


1513.0 


1813.0 


0.020 (0.314) 


0.114 (0.225) 




1 


< 7.60 


NGC5864 


3.71 


5.07 


1724.0 


2024.0 


0.070 ( 0.354) 


-0.468 ( 0.359) 




1 


< 7.74 


NGC5866 


















8.47 (0.01) 


NGC5869 


3.67 


4.25 


1915.0 


2215.0 


-0.357 ( 0.376) 


0.066(0.317) 




1 


< 7.63 


NGC6010 


3.69 


3.76 


1872.0 


2172.0 


-0.082 ( 0.352) 


-0.061 ( 0.272) 




1 


< 7.78 


NGC6014 


5.28 


5.07 


2266.0 


2570.0 


7.482 ( 0.430) 


10.732 (0.413) 


1.43 (0.10) 


1 


8.77 (0.02) 


NGC6017 


4.00 


3.66 


1638.0 


1938.0 


0.885 ( 0.347) 


0.613 ( 0.283) 




1 


< 7.73 


NGC6149 


3.80 


4.22 


2277.0 


2577.0 


0.614 (0.312) 


-0.082 ( 0.354) 




1 


< 7.90 


NGC6278 


3.34 


3.26 


2682.0 


2982.0 


-0.233 ( 0.287) 


0.325 ( 0.335) 




1 


< 7.98 


NGC6547 


3.45 


5.21 


2527.0 


2827.0 


-0.400 ( 0.332) 


0.796 ( 0.418) 




1 


< 8.00 


NGC6548 


















< 7.58 


NGC6703 


3.70 


7.12 


2223.0 


2523.0 


-0.384 ( 0.344) 


-0.158 (0.563) 




1 


< 7.62 


NGC6798 


2.61 


2.79 


2210.0 


2510.0 


0.795 ( 0.205) 


0.601 ( 0.204) 


< 0.77 


1 


7.83 (0.10) 


NGC7280 


2.91 


3.55 


1695.0 


1995.0 


0.401 (0.298) 


0.427 ( 0.273) 




1 


< 7.49 


NGC7332 


2.93 


3.44 


1047.0 


1347.0 


0.350 ( 0.283) 


0.445 ( 0.259) 




1 


< 7.41 


NGC7454 


3.31 


3.46 


1870.0 


2170.0 


-0.037 ( 0.250) 


0.247 ( 0.271) 




1 


< 7.39 


NGC7457 


3.00 


3.71 


694.0 


994.0 


-0.195 (0.299) 


0.199 (0.271) 




1 


< 6.96 


NGC7465 


5.45 


5.23 


1827.0 


2117.0 


11.887 (0.446) 


21.320(0.478) 


1.79 (0.08) 


1 


8.79 (0.02) 


NGC7693 


3.48 


4.04 


2352.0 


2652.0 


0.515 (0.320) 


0.491 (0.311) 




1 


< 7.86 


NGC7710 


3.59 


4.73 


2257.0 


2557.0 


0.070 ( 0.296) 


-0.206 (0.351) 




1 


< 7.80 


PGCO 16060 


3.28 


4.14 


2556.0 


2947.0 


2.105 (0.317) 


1.661 ( 0.363) 


0.79 (0.21) 


1 


8.26 (0.06) 


PGC028887 


3.53 


7.42 


2683.0 


2983.0 


0.609 ( 0.346) 


1.463 ( 0.642) 




1 


< 8.03 


PGC029321 


4.08 


5.76 


2747.0 


2902.0 


3.335 ( 0.284) 


4.762 ( 0.296) 


1.43 (0.15) 


1 


8.53 (0.04) 


PGC035754 


3.55 


5.17 


2384.0 


2684.0 


-0.285 ( 0.285) 


-0.353 ( 0.358) 




1 


< 7.90 


PGC042549 


3.76 


3.89 


2672.0 


2972.0 


0.203 ( 0.388) 


0.341 ( 0.295) 






< 8.07 


PGC044433 


3.53 


4.52 


2525.0 


2825.0 


0.230 ( 0.328) 


0.165 (0.329) 






< 7.98 


PGC050395 


3.21 


3.44 


2172.0 


2472.0 


0.188 (0.294) 


0.042 ( 0.274) 






< 7.87 


PGC051753 


2.79 


4.12 


2268.0 


2568.0 


0.253 ( 0.308) 


0.369 (0.314) 






< 7.92 


PGC054452 


3.60 


4.66 


1768.0 


2068.0 


-0.410 ( 0.335) 


0.088 ( 0.368) 






< 7.73 


PGC056772 


2.13 


2.69 


2508.0 


2821.0 


1.620(0.201) 


2.049 ( 0.203) 


1.26 (0.20) 




8.19(0.05) 


PGC058114 


4.52 


3.60 


1383.0 


1768.0 


11.425 (0.450) 


21.487 (0.315) 


1.88 (0.08) 




8.60 (0.02) 
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Table 4 (cont'd) 



Galaxy 


rms (1-0) 


rms (2-1) 


Vel. 


range 


-f(l-O) 


-7(2-1) 


Ratio 


Source logM(H 2 ) 




mK 


mK 


km s 


Kkms" 1 


Kkms -1 


(2-l)/(l-0) 


M 


PGC061468 


2.83 


3.31 


2282.0 


2496.0 


1.260(0.211) 


0.949(0.214) 


0.75 (0.21) 


1 8.00 (0.07) 


PGC071531 


2.39 


2.63 


1880.0 


2180.0 


0.333 ( 0.262) 


-0.116(0.206) 




1 < 7.65 


PGC170172 


3.97 


4.53 


2412.0 


2712.0 


-0.066 ( 0.373) 


0.240 ( 0.376) 




1 < 7.97 


UGC03960 


2.69 


6.08 


2105.0 


2405.0 


-0.310(0.318) 


0.031 (0.481) 




1 < 7.81 


UGC04551 


3.26 


5.28 


1578.0 


1878.0 


-0.060 ( 0.292) 


-0.317(0.501) 




1 < 7.62 


UGC05408 


2.96 


3.39 


2868.0 


3123.0 


1.633 (0.259) 


3.409 ( 0.241) 


2.09 (0.36) 


1 8.32 (0.06) 


UGC06062 


4.02 


4.72 


2484.0 


2784.0 


-0.038 (0.314) 


-0.486 ( 0.366) 




1 < 7.93 


UGC06176 


3.71 


4.79 


2540.0 


2870.0 


3.918 ( 0.380) 


8.503 ( 0.455) 


2.17 (0.24) 


1 8.58 (0.04) 


UGC08876 


3.02 


4.12 


1935.0 


2235.0 


0.450 ( 0.299) 


0.059 ( 0.306) 




1 < 7.80 


UGC09519 


4.73 


4.01 


1501.0 


1818.0 


12.662 (0.436) 


13.857 ( 0.327) 


1.09 (0.05) 


1 8.77 (0.01) 








Observed, but not members of ATLAS 3D 






IC2099 


2.92 


4.71 


2649.0 


2944.0 


0.813 (0.266) 


1.861 (0.354) 


2.29 (0.87) 


1 7.82(0.12) 


IC3102 


3.67 


8.13 


2085.0 


2385.0 


0.256 ( 0.338) 


1.736(0.613) 




1 < 7.23 


NGC2697 


4.20 


10.85 


1670.0 


1979.0 


4.365 ( 0.461) 


5.144(0.965) 


1.18(0.25) 


1 8.25 (0.04) 


NGC4292 


3.50 


7.21 


2151.0 


2397.0 


2.923 (0.315) 


2.607 (0.512) 


0.89 (0.20) 


1 7.66 (0.04) 


NGC4309 


4.79 


8.25 


951.0 


1173.0 


6.828 ( 0.390) 


7.342 ( 0.564) 


1.08 (0.10) 


1 8.05 (0.02) 


NGC4352 


2.93 


4.83 


1932.0 


2232.0 


0.048 ( 0.330) 


-0.416 (0.397) 




1 < 7.03 


NGC4479 


3.92 


5.35 


726.0 


1026.0 


-0.453 ( 0.438) 


-0.656 ( 0.391) 




1 < 7.33 


NGC4620 


4.45 


7.92 


994.0 


1294.0 


-0.277 ( 0.420) 


-0.680 ( 0.706) 




1 < 7.52 


NGC7743 


5.72 


4.80 


1561.0 


1805.0 


5.673 (0.381) 


10.746 ( 0.339) 


1.89 (0.14) 


1 8.14(0.03) 


PGC2800919 


3.02 


4.03 


2773.0 


3086.0 


1.230(0.316) 


2.495 ( 0.282) 


2.03 (0.57) 


1 8.10(0.10) 


UGC01756 


4.09 


6.02 


2862.0 


3162.0 


0.860(0.331) 


-0.536 ( 0.507) 




1 < 8.00 


UGC05467 


4.13 


6.40 


2815.0 


2986.0 


4.055 ( 0.290) 


3.399 ( 0.385) 


0.84 (0.11) 


1 8.61 (0.03) 



Note. — Intensities are quoted as main beam brightness temperatures T m b- The rms noise values are for channels binned to 31 km s~ . 
The velocity range is either 300 km s _1 (for nondetections) or the actual range used in the sums (for detections). All quoted uncertainties 
are la; H2 mass upper limits are 3 a. Where the line ratio i s a limit, it is computed using the 3cr upper limit on the undetected line. H2 
masses assume the distances listed in Cappellari et al. (201 la). For galaxies taken from the literature, we give onl y the H 2 mass correct ed to 
our assumed d i stance and other quant i ties m ay be found in the original sour ces. Da t a sour ces: 1 = this paper, 2 =ICombes et alj d2007h . 3 = 
IWelch & Sage] J2003h . 4 = ISage et alj J2007l) . 5 = IWiklind et al] i ll 9951) . 6 = lYound J2002h . 7 = I Welch. Sage. & Yound feOld) . The table is 
available in its entirety from the project website http://purl.com/atlas3d 
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Table 5. Fitted properties of detected lines 



Galaxy 


7 (i-o) 


J (2-l) 


V (l-0) 


V (2-l) 


AV(i_o) 


AV (2 _i) 


Tmb (1-0) 


?mb (2-1) 


Comme 




Kkms" 1 


Kkms" 1 


km s 


km s 


km s 


km s 


mK 


mK 




IC0676 


12.11 (0.40) 


17.03 ( 0.33) 


1429.0 ( 1.6) 


1428.0 ( 0.8) 


97.4 ( 3.4) 


87.1 ( 1.7) 


116.8 


184.3 




IC0719 


3.77 ( 0.37) 


6.10(0.33) 


1859.0(14.8) 


1836.0(6.8) 


292.6 (29.0) 


267.3 (14.1) 


12.1 


21.5 




IC1024 


12.03 ( 0.44) 


16.24 ( 0.43) 


1466.0 ( 3.0) 


1492.0 ( 1.9) 


160.9 ( 5.9) 


146.9(3.8) 


70.3 


103.9 


a;2h? 


NGC0509 


0.49(0.15) 




2288.0 ( 5.9) 




35.4(11.6) 




13.0 






NGC0524 


3.11 (0.36) 


5.30(0.52) 


2393.0 (21.2) 


2412.0(17.1) 


341.2 (37.0) 


336.6 (34.8) 


8.6 


14.8 


2h? 


NGC1222 


18.36(0.41) 


29.41 ( 0.47) 


2410.0 ( 1.6) 


2399.0 ( 1.0) 


142.4 ( 3.3) 


132.6(2.2) 


130. 


237. 


a* 


NGC1266 


32.30 ( 0.47) 


70.54 ( 2.50) 


2161.4(0.9) 


2161.9(0.4) 


129.4(2.3) 


127.2 ( 2.0) 


234.5 


520.8 








36.92 ( 2.27) 




2151.9(4.3) 




349.9 (18.3) 




99.1 


wings 


NGC2685 


1.21 (0.23) 




819.0(9.4) 




99.4 (21.2) 




11.4 






NGC2764 


15.27 ( 0.39) 


25.44(0.81) 


2707.1 (0.6) 


2702.0 ( 3.3) 


288.8 ( 1.1) 


226.0 ( 7.0) 


56 


105.8 


2h 


NGC2824 


4.49 ( 0.29) 


6.26 ( 0.49) 


2739.0 ( 7.8) 


2724.0 ( 8.0) 


235.2 (16.4) 


219.9(17.7) 


18.0 


26.8 




NGC3182 


2.92 ( 0.37) 


3.76 ( 0.20) 


2086.0(13.2) 


2127.7(0.8) 


172.3 (27.7) 


192.5 ( 1.6) 


15.9 


27 


2h 


NGC3245 


0.94 ( 0.22) 




1403.0 (10.5) 




87.4 (21.4) 




10.1 






NGC3599 




0.92 ( 0.20) 




791.7(9.5) 




91.7 (25.5) 




9.4 




NGC3619 


4.61 ( 0.45) 


3.97(0.61) 


1576.0 (20.6) 


1534.0 (21.7) 


390.8 (40.0) 


303.8 (50.2) 


11.1 


12.3 




NGC3626 


7.47 ( 0.55) 


11.51 (0.56) 


1448.0(12.1) 


1428.0 ( 6.7) 


310.3 (23.9) 


295.7 (15.4) 


22.6 


36.6 


a 


NGC3665 


12.31 (0.45) 


14.25 (0.91) 


2084.6 ( 0.7) 


2140 


678.1 ( 1.5) 


575 


31. 


52. 


a*; 2h 


NGC4036 




3.77 ( 0.32) 




1432.0(11.7) 




300.7 (24.9) 




11.7 




NGC4119 


4.70 ( 0.30) 


5.70(0.48) 


1663.0 ( 4.0) 


1667.0(4.6) 


119.7(7.6) 


115.0(9.2) 


36.9 


46.5 


2h? 


NGC4281 




0.56(0.17) 




2360.0 ( 5.2) 




36.7 (13.9) 




14.4 


2h? 






1.28(0.26) 




2994.0 ( 7.8) 




94.9 (24.4) 




12.7 




NGC4324 


3.14(0.33) 


2.33 ( 0.35) 


1667.0 ( 9.5) 


1712.0(17.4) 


182.8 (22.6) 


241.3 (34.0) 


16.1 


9.1 




NGC4429 


7.18 (0.36) 


7.08 ( 0.54) 


1110.0 ( 1.1) 


1110.0 ( 1.7) 


549.9 ( 2.2) 


501.7(3.5) 


16. 


24. 


2h 


NGC4435 


4.74(0.31) 


5.79 ( 1.05) 


775.8 ( 1.1) 


688.3 (28.3) 


440.4 ( 2.2) 


319.8(84.2) 


15. 


25. 


a*; 2h 


NGC4596 


0.37(0.10) 


1.91 (0.33) 


1774.0(4.8) 


1954.0(11.3) 


31.1 (7.3) 


137.1 (22.6) 


11.1 


13.0 


2h? 


NGC4596 


0.60(0.16) 




2005.0(8.1) 




50.2(18.2) 




11.3 






NGC4643 




0.53(0.12) 




1263.0(3.1) 




28.8(5.3) 




17.2 


2h? 






0.34(0.10) 




1629.0(4.4) 




26.6(5.8) 




12.0 




NGC4684 


1.73 (0.28) 


5.91 (0.54) 


1575.0(12.5) 


1579.0 ( 7.2) 


142.7 (20.4) 


163.0(16.2) 


11.4 


34.0 




NGC4694 


6.20 ( 0.35) 


6.73 ( 0.29) 


1188.0(2.2) 


1172.0 ( 1.4) 


80.2 ( 5.3) 


70.9(3.1) 


83. 


96.3 


a* 


NGC4710 


32.74 ( 0.76) 


40.44 ( 0.68) 


1133.0(2.7) 


1123.0 ( 1.6) 


231.2(5.8) 


209.5 ( 3.8) 


133.5 


181.2 




NGC4753 


11.67 (0.59) 


14.61 (0.53) 


1217.0(10.5) 


1238.0(6.5) 


413.2 (24.0) 


387.0(14.5) 


26.5 


35.5 


a 


NGC5173 


2.12(0.31) 


3.32 ( 0.26) 


2399.0(11.8) 


2412.0 ( 6.5) 


147.2 (30.9) 


174.4(14.1) 


13.6 


17.9 


a 


NGC5273 


1.40(0.24) 


3.44 ( 0.20) 


1088.0(11.7) 


1085.0(3.2) 


137.2 (31.3) 


117.9(6.8) 


9.6 


27.4 




NGC5379 


3.90(0.31) 


3.10(0.43) 


1797.0(4.1) 


1795.0(3.2) 


103.2(9.1) 


49.9(7.1) 


35.5 


58.4 




NGC5638 




0.39(0.10) 




1653.0 ( 3.0) 




24.9(7.1) 




15.0 




NGC6014 


8.11 (0.40) 


11.55 (0.39) 


2410.0 ( 4.2) 


2413.0(2.6) 


165.2 ( 8.0) 


155.7(4.7) 


51. 


77. 


a*; ring' 


NGC6798 


2.25 ( 0.36) 




2498.0 (35.0) 




437.5 (76.4) 




4.8 






NGC7454 


0.17 (0.07) 


0.41 (0.14) 


2017.0 ( 2.4) 


2046.0 ( 4.4) 


9.5 ( 6.2) 


30.8 (12.6) 


16.7 


12.7 




NGC7465 


12.53 (0.39) 


22.38 ( 0.47) 


1974.0(2.1) 


1974.0 ( 1.4) 


129.8(4.1) 


143.9(2.8) 


94. 


176. 


a*; 2h 


PGCO 16060 


2.09 ( 0.31) 


1.66(0.35) 


2751.0(19.0) 


2768.0 (19.7) 


251.3 (40.8) 


195.3 (43.7) 


7.8 


8.0 





Table 5 (cont'd) 



Galaxy 


7 (i-o) 


1(2-1) 


V(i_o) 


V(2-l) 


AV (1 _ 0) 




Tmb (1-0) 


^mb (2-1) 


Comments 




Kkms -1 


Kkms -1 


km s 1 


km s 1 


km s _1 


km s _1 


mK 


mK 




PGC029321 


3.58 ( 0.32) 


4.82 ( 0.34) 


2826.0 ( 4.4) 


2823.0 ( 2.9) 


98.7 ( 9.4) 


89.9 ( 6.5) 


34.1 


50.4 




PGC056772 


1.70(0.22) 


2.02 ( 0.22) 


2656.0(14.8) 


2698.0 ( 9.6) 


228.6 (34.3) 


155.2 (20.2) 


7.0 


12.2 


a 


PGC058114 


11.63 (0.34) 


22.27 ( 0.27) 


1507.0 ( 2.2) 


1522.0 ( 0.9) 


147.2 ( 4.9) 


160.5 ( 1.9) 


74.2 


130.3 


a;2h 


PGC061468 


1.32(0.22) 


1.06(0.20) 


2390.0(11.8) 


2385.0 ( 9.4) 


136.0 (21.7) 


100.5 (13.1) 


9.1 


9.9 




UGC05408 


1.48(0.25) 


3.34 ( 0.24) 


2963.0 ( 9.6) 


2972.0 ( 3.7) 


107.2 (22.8) 


103.8(8.4) 


13.0 


30.3 


a; 2h? 


UGC06176 


4.04 ( 0.37) 


9.13 (0.52) 


2713.0(10.1) 


2692.0 ( 6.4) 


207.7 (18.6) 


236.2 (12.4) 


18.3 


36.3 


2h 


UGC09519 


13.17(0.42) 


13.75 (0.33) 


1665.0 ( 2.9) 


1686.0 ( 1.7) 


179.5 (6.3) 


155.0(4.1) 


73. 


83.3 


a* 








Observed, but not members of the ATLAS 3D 


sample 








IC2099 


0.97(0.18) 


1.90(0.29) 


2758.0 ( 9.4) 


2783.0(7.2) 


94.6(16.5) 


100.8 (15.1) 


9.6 


17.7 




NGC2697 


4.74 ( 0.53) 




1784.0 (14.6) 




254.8 (29.3) 




17.5 






NGC4292 


3.10(0.33) 




2280.0 ( 9.3) 




166.8 (18.4) 




17.4 






NGC4309 


7.07 ( 0.37) 


7.92(0.57) 


1055.0 ( 2.9) 


1056.0 ( 3.8) 


110.5 (6.3) 


112.1 (7.9) 


60.1 


66.4 




NGC7743 


5.83 (0.33) 


11.06(0.34) 


1661.0(3.0) 


1661.0 ( 1.6) 


104.7 ( 6.7) 


115.6(3.8) 


52.3 


89.9 




PGC2800919 


1.30(0.25) 


2.62 ( 0.27) 


2920.0(16.7) 


2937.0 ( 8.5) 


153.1 (22.5) 


161.7(15.4) 


8.0 


15.2 


a 


UGC05467 


4.29(0.31) 


3.72 ( 0.39) 


2888.0(4.1) 


2878.0 ( 4.4) 


108.0 ( 8.0) 


84.8 ( 8.2) 


37.3 


41.3 





Note. — Columns: (1) Source name; (2)-(3) line area in CO(1-0) and CO(2-l); (4)-(5) central velocity in CO(1-0) and CO(2-l); (6)-(7) FWHM in CO(1-0) 
and CO(2-l); (8)-(9) peak brightness temperature in CO(1-0) and CO(2-l). Formal Icr uncertainties in each quantity are quoted in parentheses. Comment 'a' in 
column (10) indicates the line profile is notably asymmetric, and 'a*' indicates the line is sufficiently asymmetric that the peak brightness temperature (and area, 
systemic velocity and width, for NGC 3665 in 2-1) has been measured directly off the profile rather than with a fit. Comment '2h' indicates a double-horned line 
profile. Notes on specific galaxies: NGC 1266 - the 1-0 line is, formally speaking, adequately fit with one Gaussian component but the 2-1 line has broad, low-level 
wings which require two Gaussian components of the same center and different widths. NGC 6014 - the line profile is neither a double horn nor a Gaussian, but 
the central dip suggests a ringlike CO distribution. NGC 4281, 5638 and 7454 - the spectrum shows hints of two peaks. The peaks were fitted independently and 
both velocities are listed. Formally the Gaussian fit detects each of these with an integrated area greater than 3cr, but as there is no emission within 300 km s _1 of 
the optical velocity it is a nondetection in Tableland it is NOT counted as a detection the rest of the analysis. NGC 4643 is formally detected via the integrated 
intensity in the central 300 km s -1 of the 2-1 spectrum and its line is also suggestive of two peaks. As for NGC 4281, the two velocities listed are those of the 
two peaks (rather than the systemic velocity that one would infer from assuming those are two horns). 



